Solutions of quasi-linear wave equations 
poly homogeneous at null infinity in high dimensions 

Piotr T. Chrusciel 
Gravitational Physics 
University of Vienna 

Roger Tagne Wafo 
Departement de Mathematiques et Informatique 
Faculte des Sciences, Universite de Douala 

October 13, 2010 



Abstract 

Wc prove propagation of weighted Sobolcv regularity for solutions 
of the hypcrboloidal Cauchy problem for a class of quasi- linear symmet- 
ric hyperbolic systems, under structure conditions compatible with the 
Einstein-Maxwell equations in space-time dimensions n + 1 > 7. Simi- 
larly we prove propagation of poly homogeneity in dimensions n+1 > 9. 
As a byproduct we obtain, in those last dimensions, polyhomogeneity 
at null infinity of small data solutions of vacuum Einstein, or Einstein- 
Maxwell equations evolving out of initial data which are stationary 
outside of a ball. 
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1 Introduction 



A problem of current interest is the asymptotic behavior of solutions of hy- 
perbolic equations in the radiation zone. For large (however, not for all) sets 
of initial data, this question can be reduced to one where the initial data are 
given on a Cauchy surface that resembles a hyperboloid in Minkowski space- 
time. In recent works [12, 13], poly homogeneity of solutions of such Cauchy 
problems, with polyhomogeneous initial data, has been proved for a large 
class of semi-linear symmetric hyperbolic systems. The object of this work 
is to extend those results to quasi-linear equations satisfying certain struc- 
ture conditions which are compatible with the vacuum Einstein equations, 
or with the Einstein-Maxwell equations, in space-time dimensions n+1 > 9. 

A special case of our results is Theorem 7.2 below, where polyhomogene- 
ity at null infinity of small data global solutions of the Einstein-Maxwell 
equations, evolving out of initial data which are stationary outside of a 
compact set, is established; this is perhaps the most significant result in this 
work. For clarity we repeat the relevant part of that theorem here: 

Theorem 1.1 In dimensions n + 1 > 9 the global solutions of Einstein- 
Maxwell equations constructed in [18, 19] out from small initial data sta- 
tionary outside of a compact set are polyhomogeneous at null infinity. 

The polyhomogeneous expansions above are in terms of powers of log r 
and negative integer powers of r in odd space dimension, while one has 
powers of log r and negative half-integer powers of r in even space dimension. 

Theorem 1.1 should be compared with [4], where even space-time dimen- 
sion n + 1 > 6 is assumed, where initial data Schwarzschildian outside of a 
compact set are considered, and where solutions which are smooth at null 
infinity are obtained. The methods of that last reference completely fail in 
odd space-time dimensions. Furthermore, in odd space dimensions, generic 
initial data which are only stationary, as opposed to Schwarzschildian, are 
likely to be polyhomogeneous, but not smooth, at null infinity, and generic 
such initial data are expected to be too singular to be covered by the ap- 
proach in [4]. We also note the analysis in [3], which implies smoothness at 
null infinity of exactly stationary vacuum or electro-vacuum space-times, in 
even space-dimension, in space-time harmonic gauge. But the dimensions 
covered in [3] are precisely those not covered by the evolution theorems 
in [1,4]. 
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2 Poly homogeneity of solutions 



2.1 Notation 

The notation of [12] is used unless explicitly stated otherwise. However, 
to avoid a clash of notation with the symbol which is customarily used for 
the conformal factor arising in the rescaling of the metric, we will use the 
symbol % for the sets f2 of [12]: 

= {(x,v A ,y) : < x < y , v = (v A ) £ ff , < y < 2T) , (2.1) 

where G is a compact manifold without boundary. We will write z for the 
joint set of variables (x,y,v ). 

Let W a be a family of spaces, where a is a decay index, e.g. W a = 
*£=0},fc(*0. OT W ° = ^ < x < y}!Oo m, etc. We define 

w <a = n a<a w a . 

This notation is very useful to accommodate ln n x factors that arise in the 
problem at hand: for example, in this notation we have 

We use a slight generalization of a definition of [13]: We shall say that a 
function H(z,w) is £/^ 0<x<ry ypolyhomogeneous in z with a uniform zero of 

order I in w if the following hold: First, H is smooth in w G M. N at fixed 
z S % . Next, it is required that for all B G R and k £ N there exists a 
constant C{B) such that, for all \w\ < B and < i < min(fc, I), 

\\diH(;w)\\ o m <C(B)\w\ l -\ (2.2) 

{0<x<y},k — i v ' 

Further, 

y {0<x<y} 

at fixed constant w. Finally we demand the uniform estimate for constant 
w's: Ve > 0, M > 0, i, k £ N 3 C(e, M, i, k) V\w\ < M such that 

\\dlH(;w)\\ vr m <C(e,M,i,k). (2.4) 

The qualification "in w" in "uniform zero of order / in w v will often be omit- 
ted. Similarly to [13], the small parameter e has been introduced above to 
take into account the possible logarithmic blow-up of functions in ^o<a;<j/} 



ViGN dlH(-,w) € ^ 5 0<x<v} (2.3) 



4 



at x = 0; for the applications to the nonlinear scalar wave equation or to 
the wave map equation on Minkowski space-time, the alternative simpler re- 
quirement would actually suffice: V M > 0, i, k € N 3 C(M, i, k) \/\w\ < M 

\\dlH(;w)\yo o<x<y}km < C(M,i,k) , (2.5) 

again for constant w J s. Functions which are smooth in (w,z), and have a 
zero of order ! in w at w = 0, satisfy the above conditions. 



2.2 The theorem 

Let ip = (tpi,ip2) and set 

f:=(ip,ip), f := (ijjx,xi/)2,xip) . (2.6) 

We shall say that a function G satisfies the NL- condition if there exist 
N,pi,qi,nii E N and functions Hi with a uniform zero of order mi in the 
variable 

Wi := x qi5 (f,x 2 d x f,x 2 d y f,xd A f) 

= x q '\ipi,xip 2 ,xip, x 2 d x f, x 2 d y f, xd A f) 



such that 



G = Y,^ Pl5 Hi(z, Wl ) , (2.7) 



8=1 

with, for i = 1, . . . , N, 

V - i 

mi > £ . (2.8) 

Our first main result is the following: 

Theorem 2.1 Let be defined in (2.1), suppose that p e Z, q, 1/5 G N*, 
k £ NU {oo}, and let 

and cp, with 

& e «Sii< w },oo n *SW.o . <h,<P g «Si<»},oo > ( 2 - 9 ) 

6e a solution on % of the following system of equations: 



dyip + -B w i^ + B^tp = Ltptptp + L^tp + a + G^ 
d x ip + -B^yj + B^tp = L^cp + L^tp + b + G^, ' 



(2.10) 
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with the operators 

Lij = Lp A + xtydy + xL%d x (2.11) 

satisfying 

G xS ^{o< x <y} , G &f{o< x < y } (2.12) 

(no symmetry hypotheses are made on the matrices L^), while 

B w £ Coo(r2) + x&s ^{Q< x < y } i B^.B^.B^ G £*{o<a:<j/} ( 2 -13a) 
o,6GrW^ }) (2.13b) 

= <p e x- 1+5 ^ =0} , = ^ e x- 1+ <W£ =0} . (2.14) 

// i/ie non-linear terms Gu,, G^ satisfy the NL-condition, then 

G *{0<x<v} X ^"^{O^l/} ' 

more precisely 

^ G x 5 -V/ x=0} + x'- 1 ^^} . (2.15b) 
In particular for any r > we have 

which shows that the solution is poly homogeneous with respect to {x = 0} 
on {y > r}. 

Proof: This theorem is a generalization of the semi-linear case, Theorem 
3.7 of [12], and can be proved by following step by step the proof given there. 
A detailed exposition can be found in [25]. □ 



3 Propagation of polyhomogeneity for the Einstein- 
Maxwell equations 

Let us show that Theorem 2.1 applies to the source- free Einstein-Maxwell 
equations; we will make extensive appeal to [4]. More generally, consider a 
system of second order wave equations of the form 

pp. f fpf 



I. 



for a map / with values in M. N for some N, where r\ is the (n+l)-dimensional 
Minkowski metric. (The map / in this section should not be confused with 
the map / appearing in (2.6), compare (3.18) below.) The Einstein-Maxwell 
equations in the harmonic-Lorenz gauge can be written in this form, with 
/ := (V ~ Vfiu,^), then 

H l» . = gi iv _ = v ^ v vP(g a p - r?Q/3 ) + quadratic terms 

depends only upon — 77^, while F is a quadratic form in df with coef- 
ficients depending upon g^ u — f]av Thus, in the Einstein-Maxwell case the 
source function F has a uniform zero of order two, while the functions H^ v 
all have a uniform zero of order one. 

As in [4], and similarly 1 to [13], we use a mapping 4> : x t— > y from the of 
the future timelike cone with vertex 0, I^ x (0), of a Minkowski space-time, 
which we denote (RJ + , rj x ), into the past timelike cone with vertex of 
another Minkowski space-time, , %), defined by 

<t> : K +1 ^ ^ » y° ■= • (3-2) 

It is easy to check that 4> is a bijection from l£ x (0) onto J~ (0), with 
inverse 

<F X :y a ^x a by x a := . (3.3) 

Moreover <j> is a conformal mapping between Minkowski metrics: 

r] a pdx a dx li = n- 2 rj a/3 dy a dy^ , (3.4) 
n+1 

n := -r/^yV . (3.5) 
We work within 7~ (0) and to the future of a hypersurface 



where is a function defined on all M™ + , given by 



TiD 

where we set 



,y T0 := {y° = t }, t < , 



p = \y\ 



^(y*) 2 , x:=-\y\-y°>0, y := y° - \y\ + 1 > 
\ i=l 



In [4] one works within lt x (0), while in [13] the complement of I^ x (0) is considered. 
However, the methods of [13] apply to both situations. 
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so that 

n = -x(i- y ), d x = --(gJL_ + _) , a y = __(g|____ ; j 

and 

Furthermore the flat d'Alembertian O^y associated with the coordinates y^ 
equals 

□ -us 2 ( n ~V ( d i 9 ) | 

U n , y - 40^ 1 _ ^ _ y ^ + Cj, j + (1 _ x _ y)2 , 

where A/j is the canonical Laplacian on 5™" 1 . 

It should be kept in mind that we are interested in x small and y bounded 
away from one. 

The general relation between the wave operator on scalar functions in two 
conformal metrics transforms the left-hand-side of (3.1) into the following 
partial differential operator 

' dy a dyP v/ dx a dx^ J v v ; 

We introduce the following new set of scalar functions on 

f-n-^fo^ 1 , (3.7) 
so that the system (3.1) reads 

a0 d 2 / n-SafrmS/.ifli flfl rt d 2 / 



V 



dy a dyP 



n- n -¥ {ir*(x, f, df, ddf)^L- _ f(x, /, df)} o , 

(3.8) 

and we need to analyse the structure of the right-hand side. As calculated 
in detail in [4], if we set 

A Z := ° - ""^ - , (3-9) 

which is bounded on any bounded set of we can write 

(3.10) 



<9x^ r 11 dy a \ y a 'dy^ vvy '<9y 
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(We emphasise the occurrence of a factor of x in front of each derivative 
except d y ] that will play an important role in what follows.) Similarly 

ft 2 f eft 2 ft 

+^vxaV^y a y fi ((y - i)d y + xft x ) 2 



+ 



4VXaV»0 a yP + 2x(l - y)r? A/x ] ((y - l)ft y + xd x )}f o <jr l . (3.11) 



n — 1 a 

If we now set / o = O - ^ - /, we find: 
and 

ww = ^ l - y ^ ww ~ {n ~ l) " y)) y PaV w + w v ) 

+ ^-}l( X (l-y))^D a0 f , (3.13) 

with 

:= 2(n - 3)^7^ V + 2(*(|/ " 1))%/J • (3-14) 
Collecting all this, we conclude that 

<9 2 f f ft 2 ft 

or 1 = (^-y^i^-yf^^+^-y^y^^ + iy- 1 ^) 



ftx x ftxv v v yJJ 1 v fty x dyv v y; ' Q(AW dyM) 

+4r ?Aa r ? ^y Q /((y - + xft x f + 2(n - l)x(l - y)y Q r? a (A- £ 



l <9yM) 

+ ^nrixari^y 13 + 2x(l - y)^ ((y - + xd x ) 

+ (n - 1) [(n + 1) WA/3 y V + z(l - y)v\»] j/ • (3.15) 
The second term on the right-hand side of (3.8) is 

ff)or 1 = n^A^-(n - i)n- W)) . 

(3.16) 
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Now, A®y a = fly^, and it follows from (3.10) that the right-hand side of the 
last equation can be rewritten as 



x(l-y)) 2 x 

;i - V)) V /, (x(l - y)) V (x(y - - 2^ a y a ((y - 1)9,/ + s - (n - l)y M /)) . 

(3.17) 



As shown in [12, 13], the left-hand-side of (3.1) can be brought to the form 
needed in Theorem 2.1 by setting 

*i = f, V>2 = (d y f,d A f) , <p = dj. (3.18) 

Here dj^f = d v Af, where the u s are local coordinates on the sphere. To 
bring (3.17) to the desired form (2.7), the choice 

. n + 3 n — 3 

P20 = — r— , q 2 = — r— , 



provides the supplementary power of x needed in the arguments of F to 
satisfy the structure conditions of Theorem 2.1, provided that we choose 
1/(25) G N* in even space-dimensions; any 1/5 £ N* is admissible in odd 
ones. If we assume that F has a uniform zero of order condition (2.8) 
will now be satisfied for 



n + 1 4 A 

1712 > = 1 H < ^=^ n > 4 and ni2 > < 

n — 3 n — 3 



6, n = 4; 
4, n = 5; 
3, n = 6, 7; 
2, n > 8. 

(3.19) 

(In the Einstein-Maxwell case we have 777-2 = 2, which enforces 77 > 8.) 

Let us turn our attention to the first term at the right-hand side of 
(3.1). In what follows we will consider the following restricted class of non- 

n— 1 * 

linearities: we assume that, after replacing / by £1 2 f and changing vari- 
ables a+ — > as above, the terms H a @ takes the form 



H a/3 = G a^ n V/, !]¥+ 2 9 / ^/) , (3.20) 

with a uniform zero of order rriQ. Such a structure will clearly be obtained 
from a function in (3.1) which depends only upon /, in particular this will 
be the case for the Einstein or the Einstein-Maxwell equations, with tjiq = 1. 



10 



Using (3.15) we can write 



n+3 I n— 1 

X 2~ + ~2~ 



where F\ is linear in the second, third, and fourth argument. Assuming 
(3.20), this can be rewritten as 



n+7 

~5-Fi(fT,x 



—+ 2 f, x — +2 



n-\-7 , n — 1 -■ i n n— l|i 



dy*d y pf,X 2 +2 fy,/) 

9 !/ »^/,iV+ 1 S 1/7 /) , 



where i*2 has a uniform zero of order mi = mo + 1. With the restrictions 
on 5 as before, we will obtain the right structure by setting 

n + 7 „ n — 1 



2 ' " 2 1 

and the A^L-condition will hold provided that mi := mg + 1 satisfies 



mi > 



n + 5 



1 + 



77 



n 



1 ' 



rap > < 



' 7, 
4, 
3. 
2- 

1, 



n 
n 
n 
n 
n > 



2; 
3; 

4; 

5,6,7; 



(3.21) 



In particular the structure conditions will be satisfied by the Einstein- 
Maxwell equations in space-dimensions larger than or equal to eight. 

The hypothesis (3.20) will not be satisfied in general if H^ u in (3.1) is 
a non-linear function of / and d x nf, for then H will belong instead to the 
following class of functions (compare (3.16)) 



(3.22) 



An analysis similar to the one above shows that, for H^ ul s which are a finite 
sum of terms of the form (3.22), we will obtain the right structure by setting 



Pid 



n + 5 



qi5 



n 



2 ' ^ 2 ' 

and the iVX-condition will hold provided that mi = rriQ + 1 satisfies 



mi > 



n + 3 



1 + 



71 



71 



n > 4 and mo > < 



(7, n = 4; 
4, n = 5,; 
3, n = 6; 
2, n = 7,8,9; 
[ 1, n> 10. 

(3.23) 
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The reader should have no troubles similarly working out the conditions 
on the nonlinearity for general fTs which depend on /, d x vf and d x ^d x ^f. 
Summarizing, we have proved: 

Theorem 3.1 Let f be a solution of equation (3.1), define ipi, ip2, an d <P by 
(3.18), where f is given by (3.7). Suppose that (3.19) holds, and assume that 
either (3.20) with (3.21) hold, or (3.22) with (3.23) hold. If (2.9) and (2.U) 
hold, then the conclusions of Theorem 2. 1 apply. In particular Theorem 2. 1 
applies to the Einstein-Maxwell equations in space-time dimensions n + 1 > 
9. 

4 Towards solutions with a poly homogeneous Scri 

In order to establish existence of solutions of the vacuum Einstein equations, 
in sufficiently high dimensions, with a polyhomogeneous Scri, it remains to 
construct appropriate initial data, and show that the corresponding solutions 
are in the right function spaces. 

Recall, now, that large classes of polyhomogeneous hyperboloidal ini- 
tial data have been constructed in [2] (the emphasis in that reference is 
on n = 3 at several places, but the general results there show that the 
conformal method, starting from smooth or polyhomogeneous seed fields, 
provides polyhomogeneous solutions of the general relativistic vacuum con- 
straint equations in any dimension n > 3). There is little doubt that large 
collections of initial data so constructed provide polyhomogeneous data for 
the harmonically reduced equations of the last section, but we have not 
checked this in detail. Instead, we will follow the standard-by-now strategy 
of using initial data which are stationary outside of a compact set. So, in 
Section 4.2, we provide large classes of Corvino-Schoen type initial data with 
polyhomogeneous asymptotics on hyperboloids. One of the reasons for pro- 
ceeding this way is that small such initial data lead to global, geodesically 
complete solutions [19,20]. 

One then needs to verify that the associated solutions satisfy the space- 
time weighted regularity conditions needed in Theorem 2.1. One could hope 
that the Lindblad-Rodnianski type estimates of Loizelet [19, 20] would pro- 
vide that information. It turns out that the available estimates, for space- 
times obtained by evolving small initial data of Section 4.2, are not sufficient 
for our poly homogeneity result; this is analyzed in Section 4.3. This means 
that the desired estimates have to be derived from scratch, which will be 
done in the remainder of this paper. 
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4.1 Stationary vacuum metrics in higher dimensions 

The only way, so far, of obtaining space-times with controlled asymptotic 
behavior near i° is to use initial data sets which are stationary at large 
distances. We will outline the construction of such data in Section 4.2, but 
before doing this it is convenient to start with a short discussion of stationary 
metrics in higher dimensions; our presentation follows [3]. 

Consider a vacuum Lorentzian metric n+1 g in any space-time-dimension 
n > 3, with Killing vector X = d/dt. In the region where X is timelike 
there exist adapted coordinates in which n+1 g takes the form 

™+ig = -V 2 {dt + 9 i dx i f + g ij dx i dx j , (4.1) 
=6 " =g 

3 t V = dtO = d t g = 0. (4.2) 

The vacuum Einstein equations (with vanishing cosmological constant) read 
(see, e.g., [15]) 

\ vv*vv = \\x\ 2 g , 

< Ric(g) -V-^essgV = ^AoA , (4.3) 
k div(FA) = , 

where 

Xij = —V 2 (di9j — djOi) , (A o X)ij = Xi k Xkj ■ 
We assume that there exists a > such that 

gij - Sa = 0(r~ a ) , d k9ij = 0(r- Q " 1 ) , (4.4) 

similarly for V — 1 and 6%. A redefinition t — > t + tp, introduces a gauge 
transformation 

9 -»• 6 + dip , 

and one can exploit this freedom to impose restrictions on 9. For our pur- 
poses it is convenient to impose the harmonic gauge, \3t = 0, which reads 

diiy/d^V^Qj) = . (4.5) 

Equation (4.5) can always be achieved by replacing 9 by 9 + dvjj, and solving 
the resulting linear equation for ip, cf., e.g., [5,?] for the relevant isomor- 
phism theorems.) One can then introduce new coordinates [?] which are 
harmonic for g. 
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In space-harmonic coordinates, and in the gauge (4.5), the system (4.3) 
is elliptic, and standard considerations show that the functions gij, V and 
9{ have a polyhomogeneous expansion in terms of log r and inverse powers 
of r. Furthermore, n+l g^ u is Schwarzschild in the leading order, and there 
exist constants a*- such that 



It is of interest to enquire whether or not the logarithmic powers are es- 
sential in the polyhomogeneous expansion. It has long been know in space- 
dimension three that, for metrics which are stationary and vacuum in the 
asymptotic region, coordinate systems exist where no log r terms arise when- 
ever the ADM mass is non-zero [22]. The same property is true for static 
solutions with non-zero ADM mass in space-dimension four [3] . Now, in the 
evolution theorems used below we need all coordinates to satisfy the wave 
equation, 



and the transition from the coordinates used in [3] to the coordinates satis- 
fying (4.6) might introduce log terms: This is exactly what happens for the 
Schwarzschild metric in n = 4, which does have a logarithmic term in its 
asymptotic expansion in a natural choice of wave coordinates [4] , but this is 
the only dimension where this happens for Schwarzschild. 

In general, (4.6) is achieved by changing space-coordinates x l — > x % + 
^(x 3 ) (recall that t is already harmonic), thus solving a linear equation for 
ip l ; by standard results (see, e.g., [6]) the ip^s will have a full asymptotic 
expansion in terms of powers of lnr and inverse powers of r, and so will the 
space-time metric in the new coordinate system, when transformed from the 
space-harmonic ones. In view of the calculations in [4], this implies the ex- 
istence of polyhomogeneous asymptotics of the initial data on hyperboloids 
at as needed in Theorem 2.1. 

Rather surprisingly, in even space-dimensions larger than or equal to 
six the space-coordinates used in [3] satisfy (4.6), and so does the time 
coordinate. It follows that the analysis of stationary solutions in [3] directly 
provides wave coordinates in which no log terms occur in those dimensions. 

4.2 Corvino-Schoen data in higher dimensions 

So far we have considered metrics which are exactly stationary. Now, there 
exists a construction due to Corvino and Schoen [?, ?] (see also [9, 10], and 




Uxf 1 = 



(4.6) 
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also the more recent Reference [8], where the construction is carried out 
under considerably weaker asymptotic conditions) which allows one to glue 
exactly stationary ends to asymptotically Euclidean initial data sets. Some 
details of this construction have been presented in those references in di- 
mension three only, but the construction generalises to any dimension, as 
follows: Recall that the construction requires a family of stationary reference 
metrics which cover the whole range of asymptotic charges. In dimension 
3 + 1 this is provided by the family of metrics obtained by boosting and 
translating the Kerr metrics. In higher dimensions one such family can be 
obtained by boosting and translating the Myers-Perry metrics [21]. Note 
that the question, whether or not the reference solutions have naked singu- 
larities is irrelevant for the problem at hand because here one only needs the 
solutions at large distances. (Similarly to the Kerr family, all the metrics in 
the family so obtained have a timelike ADM momentum, and therefore can 
only be glued to asymptotically flat initial data which also have this prop- 
erty; this is no restriction for well behaved initial data sets which are spin, 
or for space-dimensions up to seven, and is expected not to be a restriction 
for well behaved initial data sets in general, but this has not been proved at 
the moment of writing of this work.) 

So let R x , ek be positive constants and consider the collection, say Cr^^ 
of general relativistic electro-vacuum initial data sets (B. n ,g,K) which are 
stationary outside a coordinate ball B(R X ) and with weighted Sobolev norm 
controlling ^-derivatives of the metric smaller than e^. Here k should be 
sufficiently large as in [4,20], and the norm should be the one described 
in those references. From what has been said this collection is non-empty, 
and contains an open set (in the topology associated to the norm) around 
Minkowski space-time. 

Now, for the Schwarzschild metric in dimension n + 1 with n > 4, and in 
harmonic coordinates, the boundary of the domain of influence of a ball is 
sandwiched between two hypersurfaces t — r = const [4, Section 5.3]. This 
remains true for stationary electro- vacuum metrics because the leading order 
behaviour of the metric coincides with the Schwarzschild one (compare [7, 
Appendix A]). This implies that the maximal globally hyperbolic develop- 
ment of all initial data in Cn Xytk contains hyperboloidal hypersurfaces, the 
asymptotic region of which is contained in that part of the space-time where 
the metric is stationary. So our considerations of the previous section apply 
to this region, leading to poly homogeneous initial data on such hypersur- 
faces. Since the leading order deviation of the metric from the flat one is 

A 71 — 1 

Schwarzschildian, the tensor field h := Q. ~(g — rj), that plays a key role 
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in our analysis, is 0(x^ n 2 ) ( n 1 ^ 2 ) = 0{x^ n 3 ^ 2 ), and in fact 

h e x^ n - 3 ^ 2 {£/ { 5 x=0} n L°°) , d y0 h e z (n - 5)/ V { i =0} n L°°) , (4.7) 

with 5 = 1 on any hyperboloid whose asymptotic part is contained in the 
stationary region. 

4.3 Lindblad-Rodnianski-Loizelet metrics near 

In this section we analyze how the asymptotic behavior of the small-data 
space-times constructed in [19] (compare [16, 17]) relates to the differentia- 
bility conditions needed in Theorem 2.1. We find that sharper decay rates 
along outgoing null geodesies would be needed for a direct proof of poly- 
homogeneity using our approach. The estimates established here are then 
combined with the results of our analysis in subsequent sections to provide 
a rather more involved proof of polyhomogeneity. 

We start by recalling some notation of [16,17,19]. Let Z denote the 
following set of vectors on Minkowski space-time: 

d a = a = 0,1,... ,n; 
Z aj 3 = x a d/3 - xpda, a, P = 0, 1, . . . , n; 

n n 

Z = '^2 X a 3 a = td t + = *^ + r ®r ■ 

a=0 i=l 

Here, as usual, xq = — x° = —t, xi = x l for i = 1 . . . , n. Let the spherical 
coordinates (r, 9 A ) be defined as 

t = x° , 

r={t^) 2 ) 1/2 , (4.8) 

i=l 

x % = roj l (9 A ), i = 1, . . . , n, 

where 6 A denotes any local coordinates on the sphere S n_1 . The vector 
fields 

L = d t + d r = d t + uj % di , L = d t — d T = d t — oj l di . 

are tangent, respectively transverse, to the light cones t — r = const. We 
note 

Zq = td t + rd r ■ 
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Furthermore, the Zij's, i,j = 1 . . . , n are tangent to the spheres S*" --1 C M. n , 
and can be purely expressed in terms of the 9 A, s. 

Let T > 0, set = (T, 0, . . . , 0), in this section it is more convenient to 
consider instead the following variation of (3.2): 

y t * = - / 7FTT x tx + T tl = -^— . (4.9) 

This provides a conformal transformation from the future causal cone cen- 
tred at in the Minkowski space-time with coordinates x M to the past 
causal cone of the origin in the Minkowski space-times with coordinates y^, 
and with conformal factor = y a y a = _^ t+ ^i +r i ■ 
To make contact with Section 2 we set 



x = -y° - p, y = y° - p + 1 where p = ( ^(y 1 ) 2 ) 1/2 

i=l 

so that 



y° = \{y - x - i) 

p=i(-y-x + l) . (4.10) 

V % = \{-V - x + l)u l (v A ), i = l,...,n 

Here w 1 is a unit vector, and the v A, s denote local coordinates on S n ~ 1 in 
the y-coordinates. One can take u/(0 ) = ), i = l,...,n; we will 

make this choice, and simply write to 1 in both x^ and y^ coordinates. 
Letting ffl s be the following family of hyperboloids, 

= {x° - s = Vs 2 + r 2 } , s > , 

we will have 

<K^) = {y° = -f} 

2s 

in particular <$>{3%\) = {y° = 

The methods of Section 2 involve the vector fields 

d 

xd x , yd y , d A = -g^j, A = 1, . . . n - 1. 

By straightforward calculations one finds, keeping in mind that p = ^ t+T yi_ r -2 
for t + T > r, 

X ~ t + T + r' ~ V ~ t + T - r T ~ 2x~2(l - y) ' ~ 2^ + 2(T 
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xd x = -\{t + T + r){d t + d T ) , 

(l-y)d y = \(t + T-r){d t -d r ) , (4.11) 
&a = linear combinations of i, j = 1, . . . , n . 

The coefficients in the equation for &a above depend only upon the angular 
variables, and a finite number of coordinate patches v A can be chosen so that 
in each of those patches the coefficients are uniformly bounded together with 
derivatives of any order. 
This leads us to 

Proposition 4.1 Let T,T > 0, t > and suppose that 

l-T<t-r<T <=> <J/<1-— |— . (4.12) 

For all k £ N, V (i, <E N x N x N™ -1 satisfying i + j + \j\ < k, and for 
any function f E C fc we /iaue 

[xSJ*^/ = J] Hp(0,y)Z I f (4.13) 
|/|<fc, zez 

^|^ 7 (0,y)|<C(i,j,/,r,To) . 

Proof: Using (4.11) one can rewrite xd x and d y as 

xd x = ~{Z Q - JZoi + T(d t + w*^)) , (4.14) 

d v = on 1 S ( Z ° + ^ + T{d t - . (4.15) 

2(1 - y) * „ ' 

V « ' =:Z 

=-fi{y) 

It is thus clear that xd x , and any of its powers, have the right structure. 
Next, the factor tpi(y) appearing in (4.15) is bounded on any compact subin- 
terval of [0, 1) (note that y = 1 corresponds to the tip of the past causal cone 
centred at the origin of the y^-coordinates). One easily finds by induction 
that 

i=i 

where the functions <pi are bounded on compact subsets of [0, 1), whence the 
result. □ 
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We wish to obtain the asymptotic behavior of the fields occurring in 
Theorem 2.1 for the global solutions 

/ := {h^,vi Ap) 

of the Einstein-Maxwell equations constructed in [19]. In order to apply 
Theorem 2.1 we need 

^ = / € V<*<,},o; fa = <&f> d *f)> * = d *f) G *SS<y},co . 



where 
Now, 



f = n-—f. 



Q = -x(l-y) 
which implies that for any oelwe have 



(xd x y(n a f) = n a ^2c( a ,i,j)(xd x yf . 

3=0 



(4.16) 



Similarly, 



(yd y y(n a f) = n a Y, c'(a, i, j, x, y )( y d y yf , d y (n a f) = n a ^2 c "(<*> i, vW 

3=0 3=0 

(4.17) 

where the functions C and C" are bounded for x in, say, [0, xo], and for y 
bounded away from 1. 

The solutions constructed in [19] satisfy the following: there exists < 
5 < 1/4 such that for t > and \t — r\ < C±, and for all / there exists a 
constant C, depending upon / and C±, such that 



\Z ] f(t,x l )\ < C{l + t + r 
{dZtffax'y < C(l + t + r 



(4.18) 
(4.19) 



where 



Now, 



8 e [dt + dr, r- l d A ) = { - 2x 2 d x , 2 ^ x ^ d A ) . (4.20) 



l+t+r l + T 

1 + — — G 



t + T + r 



i + T + r 



l + T 



1 — x — y 



for T > 0, t > 
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(4.21) 



so (4.18)-(4.19) imply 

|Z 7 /(^)| < Cx^- 5 , \dZ T f(t,x l )\ < Cx^ 4 " 1-5 
From (4.16)-(4.17) and Proposition 4.1 we obtain 

{xd^&yWjix, y, v) 
= [xd^dl&l^f^x 1 ) 

= ^ E E c(i,j,m,e,n,x,y)[xd x ] m d e y d2f(t,x l ) 

0<m<i 0<£<j 

= n ^ E E c(i,j,m,£,n,x,y) £ ■ 



Using the first inequality in (4.21) we conclude that for any < e < 1 and 
for < y < 1 — ewe have 



[xd^lydyYd^fix^y^v) < [xd x ydid2f(x,y,v) 



< Cx 



while it should be clear from (4.20) that the second inequality in (4.21) does 
not provide any new information in the coordinate ranges assumed above. 
In any case the property 

(Vl = /, V>2 = (dyf, d A f)) G V<z<,},0 » * = 9 *f G ^{O^i},- ( 4 - 22 ) 

immediately follows. Unfortunately, to apply Theorem 2.1 one would need 
5 to be an arbitrary positive number, while in (4.22) 5 is a small number 
determined by the initial data. So, as already pointed out, we need to 
derive the necessary estimates by different methods. This is the purpose of 
the sections that follow. 



5 Weighted energy estimates near a null boundary 

Let (^#, g) be an (n + l)-dimensional space-time. We consider systems of 
quasi- linear of nonlinear wave equations, with diagonal principal part of the 
form 

n s u = F(. , u, du) , (5.1) 

on a neighborhood of a null hypersurface of We suppose that the back- 
ground metric Q is a smooth function of the coordinates, of the unknown 
vector valued function u, as well as its first order derivatives. 

All calculations below will be done for a real valued function u, the result 
for a vector valued function is obtained by summing over the components. 
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5.1 The hypotheses, and the geometry of the problem 
5.1.1 The hypotheses 

We will consider the Cauchy problem associated to equation (5.1), the initial 
data will be given on a hypersurface S%. We will evolve these initial data 
to obtain a solution of our problem in a past one-sided neighborhood of a 
null hypersurface 

jY = { x = 0} 

forming the boundary, or a subset thereof, of the domain of dependence of 
o5"o- Here, and throughout, x stands for a positive function such that dx 
has no zeros on {x = 0}. We will be working in a neighborhood of {x = 0}, 
chosen so that x is a coordinate there, of the form 

V = [r , n [ x ]0,x [ x ff, 

where [to, t\ [ corresponds to the time interval, ] 0, xq [ the range of the 
variable x, and & is an (n — l)-dimensional compact submanifold of j$ 
without boundary. The coordinates will be denoted by (r, x,v), with v = 
(u" 4 )^^ the coordinates on G. We assume that d T is timelike, and we 
choose the time-orientation on j$ such that the vector d T is everywhere 
future directed. 

One can think of the set % of (2.1) as a subset of the coordinate patch 
above, compare Figure 6.2, page 78. 

On the components of the metric g with respect to the coordinates 
(t,x,v), we assume the following: 

1. We suppose that 

3e > , such that - g TT > e (5.2) 

everywhere on ~f . 

2. The components g rr and q tx can be written as 

q tt = -l + xt)°(T,x,v A ) and q tt + q tx = xI) 1 (t,x,v a ) (5.3) 
where the functions f)° and t) 1 are bounded on bounded sets. 

3. On the components g xA and g xx we assume that 

g xA =0{x) and + 2q tx + q xx = \ + 0{x) (5.4) 

and we set q xA = x\) A and g TT + 2q tx + g xx = 1 + xt), where f) and \) A 
are bounded functions on bounded sets. We further suppose that 

q tt + 2g rx + g xx > o 
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4. The vector field 

Y V 3 V := 8 T - 8 X (5.5) 

is assumed to be everywhere timelike on "V and future directed. This 
vector will be used to contract the energy momentum tensor. 

The set of functions (f) , f)^) will be denoted by f)' and 0* will denote the 
inverse matrix of the matrix (q^)- 

Remark 5.1 It follows from the above that the vector Vx (where V is the 
covariant derivative compatible with the metric g) can be decomposed as 

Vx = w (1) +/?(x)w (2) (5.6) 

where w« is causal future directed, and that there exists a constant Co such 
that 

\P(x)\<C x, \J 2) \ <C \tf\ . (5.7) 

Example 5.2 As an example, consider a conformally rescaled asymptot- 
ically flat solution of asymptotically vacuum Einstein equations in Bondi 
coordinates near Scri [23], with the metric taking the form 

Qb = e 2l3 dx®dy + xdy®dy + 2>y®dy + fj, , (5.8) 

for some functions (3 and x, and a one-form field 7. (Here y corresponds to 
the Bondi retarded time u, and x = 1 /2r is half the inverse of the luminosity 
distance r. E.g., for the Minkowski metric in any dimensions, /3 = x = = 
7.) In 3 + 1 dimensions, for smoothly compactifiable metrics, the Einstein 
equations imply, for matter fields decaying sufficiently fast, that j3 = 0{x 2 ) 
as well as 

X = 0(x 2 ), lA = 0(x 2 ), (5.9) 

with derivatives behaving in the obvious way. Equation (5.9) remains valid 
for asymptotically vacuum metrics which, after conformal rescaling, are 
polyhomogeneous and C 1 (see [14, Section 6] or [11, Appendix C.1.2]), while 
for general &^f x=0 \ H L°°-polyhomogeneous asymptotically vacuum metrics 
one has [14, Equations (2.15)-(2.19) with H = X a = 0] the asymptotic 
behaviors (3 = 0(x 2 ln N x) and 

x = 0{x 2 ), j A = 0(x 2 ln N x) , (5.10) 

for some N. Here "asymptotically vacuum" requires, for polyhomogeneous 
metrics, that the components of the energy-momentum tensor in asymptot- 
ically Minkowskian coordinates satisfy (see [14, end of Section 2]) 

Tfxv = o{r~ 2 ) . (5.11) 
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We have 

det g = — — det \i , 

which, for a Lorentzian metric, shows that fj, must be a non-degenerate 
(n — 1) x (n — 1) tensor field. It is simple to check that the inverse metric 
gtt = g a Te? a (g)tL is given by the formula 

= 49 x ®(^ + (-x+|7i;)9x-7 tl )+^, (5-12) 

with /J = [i AB 8a®9b, where //" 4S is the matrix inverse to /j,ab, 7 = 
H- AB lAdB, \"f\ 2 a = ^"(7,7) = H- AB JAJb, and <g> denotes the symmetric tensor 
product. We note 

g(Vy,Vy) = fl w = 0, 

which makes clear the null character of the level sets of y, and implies, by a 
well-known argument, that the integral curves of 

Vy = g a7 <9 a yd 7 = 5 ^d 7 = 2d x 

are null geodesies. 

Consider a new coordinate system (x,v ,t), where 

^) , (5-13) 

9, = ^ T . (5.14) 
Thus 

s" = 4(- x + |-7l^)(^ - 1^)^^ - + 4(^. - - 4^»^>(^. - + M » , 

giving 

0- = 4(- X + | 7 |2) , - = i_2(- x + | 7 |2), 0^ = -2^ B 7B ,(5.15) 

q tA = h ab ib, rr = -i + (-x + l7p, g AB = ^ B . (5.16) 

This, together with (5.10), shows that (5.3)-(5.4) hold for such metrics. 



(x,y) — > (x,t = 

so that 

Q v d - I>5 

u x 7 u x 2 r ' 
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jV = {x = 0} 



Figure 5.1: The sets ^A.r (shaded) and (the outermost trapezium). 
In this picture (but not in our hypotheses) the light-cones have forty-five 
degrees slopes, as in Minkowski space-time. 



5.1.2 The slices 

In this section we describe the sets within which we obtain our estimates, 
see Figure 5.1. Let t £ [tq, 0[ run over the range of the time coordinate r of 
the previous section. 

• Let A G [0, 1] parameterize a family of spacelike hypersurfaces S\, 
which approach {x = 0} when A approaches zero, of the form 



where o~\ is a C function such that: 

— cto(t) = i.e. So = {x = 0} 

— S\ is everywhere spacelike. 

One can legitimately raise concerns about existence of the family S\ 
with global behaviour as above when the space-time under consider- 
ation is being constructed as a solution of a Cauchy problem. While 
the aim of this work is to prove that the resulting space-time will have 
properties as in Figure 5.1, this is not known a priori. Now, one way to 
proceed is to construct the solution as the limit of solutions of linear 
equations on a sequence of metrics, each of those metrics satisfying 
controlled weighted energy estimates as proved below. In particular 
each space-time in this sequence is globally hyperbolic, with the set 
{x = 0} being part of the boundary of the domain of dependence of 
the initial surface. For each metric in the sequence a relevant family 
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S\ can be constructed using e.g. Cauchy time functions; no details 
will be given as no significant difficulties are involved. This can then 
be used to justify our estimates for each metric in the sequence, and 
for the solution. 

• By S we denote a smooth spacelike hypersurface transverse to {t = tq} 
defined by 

S = {(t,x,v a ): x = a(r)}, (5.17) 
where a is a smooth function of r such that 

< cr(n) < cr(r) < o-(to) = x . 

• H v = {(t,x,v a ); T = t, a x (x) <x< a(r)}, <^ AjT = U H v . 

ro<t<T 

• H t = Ut,x,v a ); t = t, < x < <t(t)}, W T = U H t . 

T <t<T 

Note that the boundary <9^,t °f the region ^\,t is made of four pieces, 
S\, S, H\,to an d H Ai . We recall that, for 6 G R, j S N the spaces 
S?(H AjT ), igJ(H AjT ), i^f(H A>T ) and #/(H AiT ) are defined in the appendix 
of [13].' 



5.1.3 The causality properties of the boundary 

We want to show that under the assumptions we made on certain compo- 
nents of the metric, all the hypersurfaces defined above have the nature 
which will be needed when applying the Stokes' theorem or when we will 
like to use the positivity of the stress energy momentum tensor. 

The vector Vr = V^r)^ = g^S^d^ = Q TT d T + g XT d x + 8a is normal 
to the hypersurfaces H t and H A]t , and the square of its norm is fl(Vr, Vr) = 
q tt < 0. Therefore Vr is time-like and thus these hypersurfaces are space- 
like. Their past directed unit normal is 

r, = if dp = --^=(g TT d T + d XT d x + d Ar d A ) . (5.18) 
We also note de following 

x/\sr\ V\sr\ 

that is 

rj^dx" = —^=dT . (5.19) 
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As far as the hyper surf aces S\ are concerned, the functions a\ are assumed 
to be such that the normal N = V{— x + ct\(t)} is timelike and the outward 
unit normal to this hypersurface is such that the integral of the contracted 
energy momentum tensor is negative (see (5.37)). The same remark holds 
for the hypersurface S. 

5.2 Estimates on the space derivatives of the solution 

We want to derive weighted energy inequalities for solutions of (5.1). These 
inequalities will be used to prove existence of a solution satisfying the hy- 
pothesis of the theorem of poly homogeneous solution of quasi-linear wave 
equation near scri. 

5.2.1 The stress energy momentum tensor and its properties 

The stress-energy tensor of the system (5.1) is given by 

Tf+v := V^uVyU - -Qf J ,iN cl uV a u . 

The explicit form of T °, (the component of the tensor T which in general 
determines the energy density of the system) in local coordinates system is 
given by: 

T ° = V°uV u-^V a uV a u 

= {q 00 V uVou + Q 0i ViuV u} - ^{0 OO V o nVon + 2g 0t V uV l u + Q ij V iU V jU } 
= i{fl oo (V o «) 2 -0 ii V i «V i u} = -^{-q 00 (V u) 2 + \D U \ 2 } (5.20) 

with \Du\ 2 := g^ViuVjU. 

The tensor T is symmetric and its divergence is given by 

= U 5 uV v u 

= FX7 u u when u solves (5.1) . (5.21) 

Further, one of the useful properties of the tensor T is its positivity: For 
any vectors fields v a and w a both causal future-pointing we have: 

T v »v v Wil > . (5.22) 
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Remark 5.3 In the particular frame (r, x,v A ) we will be interested with, 
let us calculate the quantity T Y := T(d T — d x ,dr) = T T T — T X T which we 
will use as energy density. Prom (5.20) we have: 

?t T = \ {Q TT {druf - Q XX (d xU f - 2Q xA d x ud A U - Q AB d A ud BU ) . 

This expression shows that in the case we are concerned with, T T T cannot 
be used to control the energy of the system near {x = 0} since the metric 
component q xx can degenerate there. On the other hand we have 

t x t = Q TT d T ud x u + Q TX {d x uf + Q rA d x ud A u , 

therefore we deduce the following expression of T Y : 

T Y = ^{q tt (d T u) 2 2 Q TT d T ud xU - (q xx + 2q tx ) {d x uf 

-2 (g xA + q tA ) 8 x ud A u - g AB d A ud B u} . (5.23) 

Now, if we set 

' A = g TT + q xx + 2q xt = 1 + O(s) > (by hypothesis) 

^A = gX A +Q Ar 



K AB 



A 

then we obtain the following decomposition of T Y 



T Y = -- 
2 




-q tt (d T u - d x uf + a ( d x u + ^ ~[" > d AU ) + { 9 AB - K AB ) d A ud B 



u 



(5.24) 

The above decomposition shows that the quantity T Y controls uniformly 
the energy of the system if and only if there exists eo > (which can be 
made to coincide with the one occurring in (5.2)) such that 

A > e , and {q ab - k ab ) ( a ( b > e Ea(Ca) 2 ; (5.25) 

the existence of such a constant follows already from our previous hypothe- 
ses. It turns out that if we have a priori bounds on the L°° norms of fj" from 
above and below, this expression can be used to control all the components 
of the stress energy tensor. In fact we have 

M = \jrd a ud„u - U/Q ap d a udpu\ < C\^\\du\ 2 < C\gl\\T T T - T X T \ ; 

(5.26) 

here the constant C depends upon eo, and is allowed to change after each 
inequality symbol in general. 
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Remark 5.4 For further purposes we note that, using the vector field d T — 
d x , the principal part of the d'Alembertian has the following form: 

9 aP d a(3 = ^(d T -d x ) 2 + 2( Q TT +B TX )(dr-d x )d x + 2B TA (d T -d x )d A 

+ (£r + 20™ + 8**) dl + 2 ( s xA + q tA ) d x d A + Q AB d A d B . 

(5.27) 

5.2.2 Estimates on the first derivatives of the solution 

We want to derive some energy inequalities for the solution u of the system 
(5.1). For this purpose, we consider the weighted energy at an instant t 
of the evolution of the system defined using the vector field d T — d x ; recall 
T Y = T t_ t 



- X 



E[u(t)} = - [ x~ 2a T ¥ —d n - l v t , x (5.28) 

JH t x 

where cT Vt,x is the measure defined on {t} x {x} x 6 by the metric g (as 
will be made precise shortly) , and a < a real parameter the range of which 
will be given later. We set 

E x [u(t)} = - [ x~ 2a T Y —d n ~ x v t , x . (5.29) 

Our strategy will be to obtain a bound of E[u(t)] from an uniform bound 
(with respect to A) of E\[u(t)]. We will apply the divergence theorem to the 
energy-momentum tensor; this holds e.g. for Cj functions u (first deriva- 
tives locally Lipschitz continuous). We want to establish the following (recall 
that eo is the constant arising in (5.2) and in (5.25), while Co is defined in 
(5.7)): 

Proposition 5.5 Let a < -\. Under hypotheses (5.2)- (5 4) and (5.25), 
there exists a constant C\ depending upon eo,Co,a such that for all 

t e [t ,ti] and u G C^ 1 

satisfying (5.1), we have 

ExHr)} < cJe x [ u ( T0 )] + £ {\\F(s)\\%, oa{nXs) + (l + H^U- + \\g%^ 

X (l + Il0lli-(H A , S ) + ll5 S Hi-(H A , s ) + II (dr ~ ^)0 8 ||!oo ( h Ais) ) ^A[n(s)]}ds| 

(5.30) 
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Proof: Stokes' theorem for the vector field A'" = x~ 2a ~ 1 T v ' x Y v on 
(compare Fig. 5.1) gives 

/ x-^T^y^dS = f V li {x- 3a - 1 T v ' i Y v }dV (5.31) 
for an arbitrary differentiable vector field Y. Here 



dV = V|detg|dr A dx A d n ~ l v , (5.32) 

where det g is the determinant of the metric g. Further, on non-characteristic 
parts of the boundary, 77^ is the unit outwards pointing conormal, and 



dS = y/\ det j\cTy , (5.33) 

with y l , i = 1, . . . , n, a system of coordinates on the corresponding bound- 
ary, and 7 the metric induced on it by the metric g; i.e. 7 = j*g , j being the 
canonical injection of the boundary into the manifold. (On characteristic 
parts of the boundary, a convenient choice of rj^ and dS will be made as 
need arises). In the case under consideration, d^\ T is made of four pieces 
Ha )To , Ha, t , together with 

S X:T ■= S x n {0 < t < t} and S T := 5 n {0 < t < r} . 

Therefore the identity (5.31) reads: 

f x~ 2a ' l T^Y v r,^dS + f x~ 2a - 1 T^Y v ^dS+ f x~ 2a - l T^Y h ' Vfi dS 

JH\,t •'H A , r[) JS\ tT 

+ I x~ 2a ~ 1 T^ l Y u r ] ^dS = I V^x^^r/TjdF . 

(5.34) 

The left-hand-side of equation (5.34) is made of four terms which will be 
labeled in their order of appearance L\, L2, L3 and L4. As mentioned 
before, we choose the vector field Y = Y^d^ to be equal to d T — d x . Once 
this choice is made, let us look at each of the terms Lj, % = 1,2, 3, 4. Recall 
that (see equation (5.19)) on Ha jT we have: 

r^cfa^ = 1 dr which implies that T^V^ = 1 {T T T - T X T } 

vls TT l Vl0 Tr l 
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and dS = y| det ^\dx A d n ~ 1 v is the surface element denoted in equations 
(5.28) and (5.23) by dx d n ~ 1 vt )X - Since rjo^det q = y/det 7 on H^ T , we 
obtain that (remember that r/^d^ is past directed) 

L\ = -E x [u(t)) . (5.35) 

From this, the sign coming from the Stokes' identity shows that 

L 2 = E x [u(t )} . (5.36) 

On the hypersurfaces S\ and S, since the unit outward normal is also 
past directed and the vector field Y"d v = d T — d x future directed, we deduce 
from the positivity of the stress energy tensor that: 

L 3 < and L 4 < . (5.37) 

We can now rewrite (5.34) as: 

-E x [u(t)} + E x [ U (t )} + L 3 + U= ( V^x-^-^/F"}^. (5.38) 
Now, let us consider the right-hand side of the above equation. We have: 

= x- 2a - x [(y^)Y v + T^V^) - (2a + l^r/H^)} 

-(2a + l)x- 1 V M x {T^-T^}} 
=: Ri + R 2 + i? 3 , (5.39) 
where 

are the Christoffel's symbols of the metric q. From (5.21), we have: 

x 2a+1 \R 1 \ = \F\\V u uY u \ = \F\\(d T u-d x u)\ < - {f 2 + {d T u - d x u) 2 }) 

< c ( eo )(F 2 + |T/-r^|) . 

(5.40) 

As far as the second term is concerned, we have: 

= ^deaw = —Tides'". 
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Thus, replacing successively in the above expression 9 with r and x and 
subtracting the two expressions we find that 

x 2a+1 R 2 = ~T/ Sva (8 T - d x ) . 

From (5.26) we obtain: 

x 2a+l \R 2 \ = \T ua (d T - d x )sT\ < (n+l)C| fl «| (|g| 2 + \{d T - d x )^\ 2 ) \T T T -T X T \ . 

(5.41) 

For the third term we have, keeping in mind (5.6): 

x 2a+1 R 3 = -(2a + l)x~ 1 T/V M xy iy 

= -(2a + l)x- 1 g^T, (T V^y y 
= -(2a + l)x~% u V fl xY u 

= -{2a + l)x~ 1 T ilu Y u J 1 ^-(2a + l)^T, u Y u J 2 ^ for a< -1/2 
> „ — ' x 

>o 

> _( 2a + i)M T ^y^(2)M = _ {2a + 1} M _ T/jgt) w (2)„ 

> -C(a,Co,n)|^|(l + | | 2 + | «| 2 )|r r --r^|. (5.42) 

Let us justify the last inequality. In other words let us show that the ex- 
pression T^ T — T^ x is controlled by \T T T — T X T \. We have: 

\T^ T - T^ x \ = \d^u (d T - d x ) u - - (g^ - g MX ) Q afS d a udpu\ 

< (d^u) 2 + [(d T - d x ) u] 2 + (|g| 2 + |g»| 2 ) [8^d a ud p u) 

< C(e )(l + \ 5 \ 2 + \ Q t\ 2 )\T T T -T x T \. See (5.24) 

Inequalities (5.40), (5.41) and (5.42) show that the right-hand side of (5.39) 
can be estimated as: 

R1+R2+R3 > -C lX ~^ { (l + |f)«| + | fl »|) (l + |g| 2 + |gf + I - a,)gf ) \T Y \ + F 2 } , 

(5.43) 

where C\ = C(a,eo,Co,n) . Now from (5.38) we have 

-E x [u(t)] + E x [u(t )} + L 3 + L 4 = R!+R 2 + R 3 , 
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thus, using (5.37), we obtain the following: 

E x [u(t)] < E x [u(r )} + CiJ^J^ x- 2Q {(l + |h s | + | s |) (l + | | 2 + |g» 



2 



+ \(d T - d x ) Q *\ 2 ) \T Y \ + F 2 (s)} ds—cT^u . 



Therefore, there exists a constant C\ > depending upon n, eo , a and Co 
such that 



ExHr)] < cJe x [u(t )] + j* {ll^WII*^,.) + (l + \\$\\l<*> + l|0 tt lk~ 

X (l + ll0ll£~(H A , s ) + ll5 S Hi-(H A , s ) + II (Sr - ^)0 S |||^ ( H A , T )) E X [u{s)]}ds 



(5.44) 



and the proof is completed. □ 



5.2.3 Estimates on the higher space derivatives of the solution 

To proceed further, we would like to have an estimate similar to (5.30) 
on space derivatives of the unknown function in equation (5.1). For this 
purpose, for k G N, (3 = (/3i, /3 2 , . . . , /3 r ) G N r , with \/3\ < k; we set: 



03) 

x 



rp i_i _— 2a 



-1+2/?! jv^uV^u - ^/V^uV^U j , 



where a < —1/2 is the real parameter of the previous section, 
X^X^ 2 ■ ■ ■ Xr r , with the Xj's being the vector fields defined in [13, page 

r 

51]: for i = 2, . . . ,r, Xi = ^ Xf(v)dji, where the Xf^'s are smooth func- 

A=2 

tions bounded on bounded set with all their derivatives, and X\ = d x . Since 
the operator V is linear, as in (5.21), we have 

03) V (r) G 8 ) 

V^T/ = x 2a - 1+2 ^U 3 {^u)V v {^u) + (-2a - 1 + 2/9i)-^=^T I / . 

Now 

n s (^u) = & p (n a u) + p , & p ]u = + p , S fi ]u , (5.45) 
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for any solution of the equation (5.1). Thus 

O 3 ) ( v (t\ 3 ) 

V M T/ = x -2«-i+2/3i + p fl) ^/?] n | V|/ (^ u )+(-2a-l+2/3i)-^r i/ ' i . 

(5.46) 

Similarly to the previous section, we set: 



08) v (0) 
T Y = TJ -TJ 



X 1 



r (pi f \p> 

E k a [u(r)}=J2 -T Y dxd n ~ l ^ x and / -T 

101=0 Ht |/3|=0 Ha ' t 

(5.47) 

08) 

Remark 5.6 From (5.24) we deduce the following decomposition for T : 

r )y = ~l-Q TT (x- a -5+fr®P(d T -d x ) U y 



Y dxd n 1 u tx . 



+ {e AB - k ab ) d A (x- a -^®Pv\ d B [x~ a -\^^v^ \ . (5.48) 



Since the coefficients of the terms arising in commutating 8a and 3>P are 
uniformly bounded, from the above we find that the energy of order k con- 
trols the J^. a -norms of the first order derivatives of the unknown function 
u. That is: 

||(^-^)«||^a^ ir) + ||a x u||5 r a (HAir) +5^ ll^||^a(H A>T ) < C E k,\[ U ( T )] ■ 

A 

(5.49) 

Let us set 

1 

V^del^ 



T" := -fl^r^ = 7=— ^ ( V^ditiifl"" ) • (5.50) 



Let us define 



M{r) := \\F\\% s{Uxr) + \\(Q,(d T -d x ) Q i)\\l x{HK 



+||(0«,b«,T)||^ (H ^. (5.51) 



We claim that: 
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Proposition 5.7 Let A > 0, k e N and suppose that a < — |. Under hy- 
potheses (5.2)-(5.4) and (5.25), there exists a function C2(eo,Co, a, k,n, M). 
monotonously increasing in M, which we write as C2(M), suc/i that for all 

T G [t ,ti] 

and for all u satisfying (5.1) we have 

E k a Mr)\ 

< E k - X [u(r )} + jTc7 2 (M( S ))|^ A [«( S )] + ||^( S )||^(H A , T ) 

+ 11 ((Sr - d x )u,d x u,d A u)\\% ?{UxT) x ||(fl lt ,f)',T)||J,o (H3kiT j|da. 

(5.52) 

Remark 5.8 The reader should note that C2 does not depend upon A. 

Proof: If the right-hand side of (5.52) is infinite there is nothing to prove. 
Otherwise, the calculations that follow should be done assuming smooth- 
ness of u, and the inequality for general u's can be obtained by a density 
argument. 

The equivalent of (5.34) for space-derivatives of the solution of (5.1) 
reads: 



Y / T^Y'^dS + £ / T u »Y» Vfl dS + Y / T^Y^dS 



+ T u »Y»^dS=Y VjT/rU (5.53) 



|/3|=0 ° e ' T \p\ 
which gives the following equation: 



k r (/?) K r (f3) 

-EgMT)]+E k y[u(To)} + Y T^Y^dS+Y T v »Y" Vfl dS 

\ r>\ ^ J S\ r \ r>\ n J Sf) r 



L3+L4 < 



k r- (/3) 

Y / V^{T^Y v }dxdv. (5.54) 



\P\=o 
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Again as in the previous section we take Y u d u = d T — d x , then the divergence 
in the right-hand side of (5.54) reads: 



, 09) , W OS) 



.T 



-2a- 



-1-2/3! + ^/3] u j ^ _ ^) 



+ T/(r^ T -r^) + (-2a + 
= : i?i + i?2 + -^3- 



' 09) 09) 



(5.55) 



If we repeat the calculations in the previous section that led to (5.41) and 
(5.42), we obtain that there exists a constant C = C(n, k, Co, a, eo) > such 
that: 



\R2\<C\gt\(\ e \ 2 + \(d T -d x ) Q t\ 2 )\T 



(rf) 



(5.56) 



and, keeping in mind that the term with the worst power of x can be dis- 
carded because of a favorable sign, 



1 + |0| 2 + |0 B | 2 )\T Y 



(5.57) 



As far as the term R\ is concerned, from the inequality ab < ^(a + b ), we 
have: 



x 



2o+l-2/3i 



\Ri 



\{^F + p g , @P]u} {d T - d x ) {®P U )\ 

(d T - d x ) (Vu) 



2 

U) + 



03) 



< (^F) 2 + C(e )|T y | + 



□ n 



(5.58) 



From inequalities (5.56)-(5.58) and the fact that L3 < and L4 < we 
obtain that: 



E k a Mr)]-E k a x [u(r )} < C 



TO 



i + l|ii tt IU~ + ll0 a IM (i + ll0lli- + ll5 8 lli 



M\&-&)bHI-)e£A<s)] + \\f(s) 



l^ fc Q (H AiS ) 



ds 



k 

+ £ 

1/31=0 



-2q-1+2/3i 



to JH A 



([□g,^]u) 2 (,s)cfe di/^ds 
(5.59) 
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with C = C(n, a, k, Cq,€q). Now, let us estimate the last term of the right- 
hand side of the above inequality. From the definition (5.50) of T v we have 

□ fl = 0^ + T^, (5.60) 

and then 

[D a ,^]u = ^[d a d^^]u-T»[^,d u ]u-[^{r»d u u)-T v ®P{d u u)} 

= : A 1 +A 2 + A 3 +A 4 . (5.61) 

To estimate the first and second terms, we use the explicit form of the 
differential operator 9 : = d^X^ 2 . . . X? r = dfcxfr. Since d T and d x 
commute with S>^ ', we have (see (5.27)) 

Ai = gr[d^d a ,^]u = 2 ^[(a r -^)aA,^ /3 ]w+2( a;A +0 TA )[a a; aA,^ /3 ]w+0 AB [aA5 B ,^]n, 

and since 

Q rA [(d T -d x )d A , ®P]u = g TA d^d A X^[(d T -d x )u]-Q TA d^X^d A [(d T -d x )u] 
we obtain that (see (5.49): 

I" a .-2a- 1 +2ft U A [{d T -d x )d A ,$> } uYdxdu < c\\Q%„p x J\(d r -d x )u\\%». 

Similarly, we have 

{Q xA + 5 tA ) [d x d A , @P]u = {q xA + 5 rA ) (d A ®P{d x u) - ^d A {d x u)) , 
which leads to: 

^ 1 +^{( B * A + f A )[d x d A ,@P]u}\s)dxdv < CH^IIic^^ll^ll^a 

< cUX^ ^Mr)} . 

Similar calculations give: 

/ x- 2a - 1+ ^(g AB [d A d B ,^]u) 2 (s)dxdu < cll^ll^^^ll^H^ 



X 
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We obtain thus the following estimate for the first term of the identity (5.61): 
I x -2.-i+2^ A 2 dxdv < C||fl*||i. (HAiT) ^ A [u(r)] . (5.62) 

Again since d T and d x commute with if we develop the second term of 
(5.61), we find that: 

A 2 = T u [@P,d v ]u = T A [^,d A ]u 
and we then have the estimates: 

(5.63) 

As far as the third term is concerned, we write 

A 3 = (T v d v u) - T V @P (d v u) = 9? (T T (d T - d x )u) - T T ^ ((d r - d x )u) 

+®P {{T x + T T )d x u) - (T x + T T )^ (d x u) 
(T A d A u) - T A ^ (d A u) 
=: I + 11 + III . 

Now we will use the weighted Moser-type inequality (A. 35) of Proposition 
A. 3 of [13] to estimate the terms of A3. Its first term gives the following 



[ x -2«-l+2/?i {J}2 dx dv 



\\^& (T T (d T - d x )u) - s*T^ ((d T - d x )u) ||^ (Ha , t) 

-t- Il2 , M /o o \„.i|2 iiirri|2 



< c s ||(5 r -a>||^||TH|io + ||(a T -a>||^ % ||T-||4o x 

\ u re ft— 1 {x— 0},1 

< c(\\(d T -d x )u\\%,\\T% + \\T% ? E k a [n(r)]) . (5.64) 

\ u k {a; — 0},1 I 



For the second term: 



[ x -2«-l+2/?i {JJ}2 dxdv 

= (T« + T T ) d x u - x^{T* + T^iMf^^ 



< c s \ \\d xU \\%J\r* + T-n^ + iia.ufc ||t* + t „ 

\ u k K — l {x — 0},1 

< cf||MI^I|T" + T T ||L, + ||T" + r-||J ? ^Kr)A . 

V u fe {a; — 0},1 / 
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The same holds for the third term of A3: 

f x~ 2a ~ 1+2h {III} 2 dxdv = \\x^@P (T A d A u)-x^T A $ p (d A u)f^ {YlxT) 

< C s (\\d A u\\% a \\T A \\ 2 +\\d A uf jea \\T A \\% ? ) 

\ u k K 1 {a; — 0},1 / 

< cf||a A u||^||T^|| 2 + ||T^|| 2 E k a [u(r)}) . 

\ k {^=o},i / 

We then obtain the following estimate for the third term of equation (5.61) 
[ x- 2a - 1+2l3l {A 3 ) 2 dxdu 
< ||(d T -d>|U*||TT|||o + \\d x u\\% S \\^ T + ^%o + ||^n||^||T A |||o 



+ IITlL, + \\r* + r T \\ 2 + \\r A \\ 2 \e% x [u(t)]. 

(5.65) 

In order to estimate the fourth term of (5.61), we need to look separately 
at each of its components as we have to make sure that every d 2 comes with 
a factor of x. We write 

A 4 = A 00 + 2A TX + 2A tA + A xx + 2A xA + A AB , (5.66) 

where the labeling A ab corresponds to the terms obtained when in A\ we re- 
place Q a ^9 2 a p with its expression as in (5.27). Now we use again the weighted 
Moser type inequality of Proposition A. 3 of [13] to estimate these terms. We 
have: 

f 3.-20,-1+2* {A AB } 2 dxdv 
= WX^^ ( Q AB d A d B u) - X^Q AB ^ (d A d BU ) \u«(a Xjr) 

< C s J2(\\dAu\\% ? \\ 5 %o + \\d A u\\ 2 U% ? ) 

' J* y 1 k K {a;— 0},1 J 

< C (^\\d A u\\% ? y\\ 2 +\\g% ? ElMr)}) , (5-67) 
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and 

f x - 2a ~ 1+2Pl {A TA } 2 dxdu 

= \\x^^ (Q rA d A (d T - d x )u) - x^ S tA @P (d A (d T - d x )u) \U o c 

< C s (||(a T -^)n|| 2 ^|| «||2+||(9 T -^) U ||^||g»||L ) 

\ 1 k K {a;— 0},1 / 

< C s (\\(d T -d x )u\\% ? \\Q% + \\g% ? E k,\[ u ( T )]) ■ (5-68) 
Continuing in this way we have: 

[ x- 2a - l+2 ^{A xA } 2 dxdv 

< \\X^®P { { Q XA + Q TA ) d A d x u) - X^ { Q XA + Q TA ) ^d A d x u\\^ 



< C \\d A d x u\\%,\\ ( Q xA + f A ) §o + \\d A d x u\\ 2 ^J ( Q xA +Q rA ) 



{i=0},1 



< C^f ||^||^||(0^ + 0^) ||L + \\d x uf^\\{Q xA +^ A )\\% f ) 

' ^ y 1 k « {x— 0},1 J 

< C^f||^||^||(0^ + 0^)||L + ||(g xA + g TA )Hij, £y[«(r)]). (5-69) 

^ \ 1^—0}, 1 J 

We recall that g rr + g XT = xr/fr, x, v A ), we then obtain the following ex- 
pression for A TX . 

= ^ [r^a^ - d x )u] - fj 1 ^ [xa x (a T - cyu] 

+ b 1 ^ [^.(a r - cy«] - xf) 1 ^ [a x (a r - d x )u] . 

V v ' 

= /3!f)i®/5(a T -&)« 

Since 

|2 
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we have 



'H A)T 

< C s (\\(d T -d x )u)\\% 



,1||2 



, 1 1 1 2 



+ \\{d T -d x )u\\*.. 

k fe {x = 0},l 



< c [\\{d T -d x )u)\\%«\\nu + W\\U eiMt)] 

\ 1 73 k B -{>=0},1 ' 

On the other hand, since g TT + 2g TX + q xx = 1 + xf), we have 



(5.70) 



+ 



< a 



^0 Q ( H A,r) 

{xfrgP (xd x [d x u]) - X^X^ {d 2 x u)} ||^ (Ha , t ) 

V v ' 

(\\xd x [d xU }\\ 2 \M 2 +\\xd x [d xU }\\ 2 HbllL x 

\ u k K — l {x— 0},1 



< a. 



%^%?\\Wy° + M^ x=ohl KM T )n ■ ( 5 - 71 ) 



We note that \\x j d x ^\\j^ a < ||$|| which can be shown by induction. In 



order to estimate the term ^4 00 , we proceed as follows: 



,4 



oo 



\q tt (d T -d x ) 



& ( [-1 + xt)°] (d T - d x f u) - [-1 + xh ]^ (a T - a. 



& ( [xf)°] (0 r - a x ) 2 - [xf> ]^ (0 r - a x ) 2 u . (5.72) 



Now using equation (5.1), (5.60) and (5.27), we obtain the following expres- 
sion of (d T — d x ) 2 u: 

(d T - d x f u = -2 (q tt + 9 TX ) (d r ~ d x ) d x - (5 TT + 20™ + D 0* - 2 ^ (d T - 3*) d A 
-2(g xA + Q TA )d x d A -g AB d A d B -T°d a u + F. (5.73) 

Here the hat means multiplication with 1/q tt (recall \g TT \ > eo > ). We 
will need the following: 



Lemma 5.9 Let 

d=(xd x ,d A ), fceN', fle 



< 



(5.74) 
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We have the following estimates: 

U\\ v e <I|H| e (5.75) 

{x — 0} ,0 £q {x — 0},0 

<rlW*f n , +4ll^ )IU°°llV'll^ olo , (5-76) 

{2=0}, 1 £q {a:=0},l {1=0}, 

and 

WtWjee < ^WWje* + C(e , ||f)°||z~)|M|^ (l + ||f)°IU fc -i) , (5.77) 
with identical estimates with ^ x=0 y replaced by SSq and J^® replaced by 

PROOF: The first inequality is obvious. Next: 

c^a < \\x~ 9 lb\\L°° + \\x" 9 d \ — !— lb 1 llroo 

< — ) + — fi^HlLoo 
e *{*=o}.o 1 TT J 



+ || 5(a;f) u)|| iooW 

£q {a: = 0},l £^ {a; = 0},0 

On the other hand, from inequality (A. 27) of [13] we have: 
1 11 

< ^IMIjtf + Wag 11-^: ll^o. (5.78) 

Now, by hypothesis we have, 

1 1 _ ^ xi)°(r, x,v A ) 



q tt (t,x,v a ) — 1 + xt)°(r, x, v A ) — 1 + xf)°(r, x, v A ) 

= -l + G(T,x,v A ,xt)°) , 

where G is any function which takes the correct values in the range of 
interest, e.g., 

G(t,x,v a , P )= with such that x (p) = ^ ^ P ~ 1 ~ ^ 



-1+p * v ' AUV \ if p>l-^ 



l 
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Recall that hypothesis (5.2) reads xt)° < 1 — eo- We have (note that the 
space of functions contains constant functions) 



1 



— Ila-o < ||l||^o + ||G(. ,xl)»)yo < C (1 + ||G(., xi) 



(5.79) 



The function G satisfies the following, for any p£R: 



\\G(.,p)\\ v o =\\G(.,p)\\^o <C(e ) 

{x — 0},k {x — OJ-,0 



and for i = 0, 1 ; 



#G(. ,p) 



< C(e ) 



1 1— * 



'{a;=0},fe-i 



These two inequalities show that G has a uniform zero of order 1 at p = 0. 
Therefore, we can apply inequality (A. 31) of [13] and obtain that 



\h v 



||G(. ,xf) u )|^o <C( eo> ||f) u |Uo 
This implies (see (5.79)) 

|| — ||^o < C(e , ||b°|U°°) (l + ||f) ||^-i) 
and (5.78) leads to (5.77). 



(5.80) 
□ 



If we insert (5.73) into equation (5.72), we obtain seven commutators 
which we label , a = 1,... ,7. These terms can be estimated in the 
same way as we did before, using (A. 34), (A. 35) of [13] and Lemma 5.9. 
They will be analyzed in the order 7 — 3 — 5 — 1 — 2 — 4 — 6. Let us estimate 
the term ^4^° containing the source term F. We have 



-2a- l+20i / /l00\ 2 



{Affdxdv = \\x^@? ([xf)°]fj -x^\x^\9^Ff^ a 



- C ( l|-F|l|g« Ik^ll^o + ll-^ll^ i Q _ 1 ll X ^ |l|' { ;i:=0 - 

< C(eo)||F||^a||xf) ||L + C(e )M°||L 

U K \X — 0},1 



x { 11*11^ + II^II^C-dl^llLoc) (l + Hl^J} . (5.81) 
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The third term can be estimated as follows: 
( x- 2a - 1+2 ^{Affdxdu 

= 2\\& (xt) Q TA d A (d T - d x )u) - Xt) ®' 3 (Q TA d A (d T - d x )u) \\\ a 

< cw A d A {d T - d x )u\\%« ||xf)°|||o + c\\Q rA d A {d T - d x )u\\ 2 ^ 

U fc fe— 1 {x — 0},1 

< Cll^lllooll^-^^II^Hxf, !!^ 

+C\\X^ ? x {\\(dr-d T ) U \\% ? \\Q% +\\(d T -d x )u\\*4Q% ? } 

{x = 0},l I 1 fe — 1 K {x — 0},0 I 

< c^l^ll^lK^-a.H^iix^ii^+c^lkf, !!^ \\sHU\(d T - d x ) u \\* 

1 fc {a;=0},l fe 

+ ll^°llio l|(5 r -rH^{||0«f o +||0 S H!ooG(l|f) O |l^)(l + ||^ O |||o )}. 

{x = 0},1 1 ^ fe— 1 \ fe— 1/ J 

(5.82) 

A similar analysis applies to A® . 

As far as the first term A® is concerned, we have 

~Af = 9? (xtftfixdjidr ~ d x )u) - Xt)°^ (t)\xd x )(d T - d x )u) . 

Using again the weighted Moser-type inequality of [13], we can evaluate the 
square of its norm as follows: 

f x -**-Wfr(Af?dxdv 

= 2^®? ([xi)%\xd x )(d T -d x ) U ]j -X^[X\?\^ (^\xd x )(d T - d x )uj \\%,a 

< C ^(xd^dr - ^)n|| 2 ^||xf,°§o + \\h\xd x ){d T - d^u^Jx^^o^ J 

< c(e ) fll^llL-ll^-a^H^ii^ii^o + ll^^^-a^ii^ Jxlj ||L, ) 

\ 1 fe fe— 1 {x — 0},1 / 
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Using now inequality (A. 34) of [13], we have (the last inequality is obtained 
by using (5.77): 

\\^(xd x )(d T - d x )u\\%, a 

fe— 1 

< c (\\xd x (d T - d x )u\\%z\\i)% +\\xd x (d T -d x )u\\%, a p 

\ u ^k-l k—1 

< c\\(d T - d x )u\\%* Hb 1 |&„ +c(e )||f) 1 ||!o O ||(a T -a>||^ c 

1 ^k~l k 

< cMHb 1 !!!- ||(a T _ d x )u\?se? + c(e )\\(d T - d x )u\\% ? 
x {lib 1 11^ + c(||f, || L ^)||i, 1 ||| oo (i + ll^ !!^)} , 

which gives 



f x- 2a - 1+2 ^{Affdxdv 



< Cieo^l^idr-d^ulll^Xo + C(eo)M°||L> x II^Hi-IK^ - d x )u\\ 2 

1 k {x— 0},1 fc 

+ c( eo )H ||^ o} J|(a r - ^uiiij. {ll^ll^ + ccii^iiL-jii^iii-a + 11x^11^)} 

< C(e )||xh°||^ ^JI^IIioclKSr-SxH^a 

{x— 0},1 re 

+C(e ) (l + lkf)°||^ =0}1 ) 11(5. - d x )u\\% ? + C(||t) ||^)||h 1 ||| oo (l + \\xi)%o)} • 

(5.83) 

Since the terms A® and have same the structure, to estimate the second 
one, we just have to replace in the estimate on A® , \\(d T — d x )u\\ 2 ^, a by 

\\d x u\\ 2 jg, a and H^Eo by \\q tA + g xA \\ 2 . 

re fe — 1 fc — 1 

We continue with the most dangerous term A® - We have (recall that 

i = i/0 TT ) 

-Ag = 9? { [Xfj°](l + Xf))dlu) - [Xt)°]^ ( [1 + z(j]fig 



/ x- 2a ~ 1+2 ^ {A° 2 } 2 = \\x^^([^][l + xt)](xd x )d x u 



-x^x^hi + xWlu) \\ 2 







< (slj (i.flgu)) -x*[x&°]^ (i.^u) 

+ ||x /,1 /3 (xf)%xd x (d x uj\ -X ft [xb°]^ (^^(Sa-ti) 



(a) + (6) . 
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Now, estimating these two expressions as we did with A® , we obtain the 
following 



(b) < c UMxdjdxuYa* H u ||£o + \\\){xd x )dM\ z ^jm^i_ o}i 

1 k re— 1 {x— 0},1 

Now equations (A. 34) of [13] and (5.77) give 

W^xd^uW 2 ^,^ < C I \\{xd x )d x u\\% s \\t)\\ 2 cgo^ + \\(xd x )d x u\\ 2 s g,c Li m 



{z=0},0 

< C(e )||f)||i=cRufc 



k 



+C(e )\\d x u\\% ? {\m\l Li + ||&||i-C(||b°||L-)(l + \\xt)% Li )} 



which gives the desired estimate for (b). 

In order to estimate the term (a) we write /3 = (fti,P v ) and ^ = 
d^JSf * , with Xl v = X^ 2 ... X? r : 

& (t)°x(i.d?u)) - [xh ]^ (l.dlu) = 9? (t)°x(i.d?u)) - f) ^ {i.xdlu) 

(i.xdlu) - [xb ]^ 3 (i.d 2 x u) 
=: (l) + (2). (5.84) 

We have 

x^(2) = /3ih°x^5f9f 1 - 1 (i.^ M ) 

as well as 

x _ 2Q _ 1+2/3l(2)2 = /32([)0)2x _ 2(a _ 1) _ 1+2(/3l _ 1) ^ } 

This identity leads to 

||^(2)||^a<C||b ||i-||i.^tt||ija-l 

II 2 / II o2„,||2 ||i||2 , ||p2 i|2 



< C||fj u ||i« Ml^ll^lllll^ + II^II^Hlll^o 

< c||b°|iicc (\\dM\%«\\H%2 + -HmO . 

Using again (5.77) we have: 

<-||l||^o +C(||h || L ^)(l + ||x 2 h°||2 

fc l £q fc — i V 
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that is 

„2e,0||2 



< C(\m\ Lao ) (1 + \\xX\\^oJ • (5.85) 
Thus, 

(5.86) 

As far as the first term of (5.84) is concerned, we have: 

11^(1)11^ = ll^ 1 ^ 3 {b°(i.(xd x )d x u)) - {i.{xd x )d x u) 

< C (||i.(x^n||i a ||f) ||L + ||i.(x^)^u||^ Q ||f,°||L i 

u k ft— 1 {x— 0},1 



< C(e )^\\d x u\\% ? 

+ x {\\{xd x )d x u\\j%«\\l\\U + \\{xd x )d x uf^ a \\l 

{x=0},l I u ft— 1 ft— i 



< CM^II^II^II^HL + II^IIL \\\d x u\\% a \\l\\l + 

1 k {x — 01,1 I 1 ft— 1 



'■O=o}, 

< C(e )\\\d x u\\% ? \fr% + \\^o E^ x [u(r) 

+^o)l|f)°ll^ o} J\d x u\\%„ {c(||f)°|M (i + H ||^oJ} 

Equations (5.86) and (5.87) finish the proof of the desired estimate for (a), 
and hence for Aq . 

Now let us consider the sixth term Aq° of A 00 . We have 

r><9„ = f T (d T - d x ) + (? x + r T ) d x + f A d A , 



and we decompose A® as 

Af = a + b + c. (5.88) 

We have 

a := ^ ([xf) ]T T (<9 T - d x )uj - [xt)°]^ (f T (d T - d x )u^j , 
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and 

x- 2a - 1+2 ^a 2 dxdv 
= \\x^@P {[x\)°}T T {d T - d x )u) - x^[x^]^ (f T (d T - d x )u) \\lg, 

< C 

< C(e )\\T r \\l x ,\\(d T - d x )u\\%J\xt)Yc?o + ||xf, L „ 

x\(d T -d x )u\\% s \\f%o + \\(d T -d x )u\\ 2 \\f T \\ 2 x 

u fe— 1 K — l {x— 0},0 

Now, from (5.75) and (5.77) we have 

HT-IU <C(eo)irr-|| £ oo, 

{x=0},0 

and 

||T%o <C(e ){\\T% +||THl!o O C7(||f, ||^)(l + ||xf, ||^ )} • 

k — l \ k — l \ k — 1/ ) 

Thus 

l x -2 a -l+2 Pla 2 dxdu < C(€ )||T T ||i oo ||(a r -5 a! )u||^||x^E 

+C(e )\\xi) \\ 2 x \\T^\\(d T -d x )u\\ 2 

\X = 0},1 K—l 

+C(||^||^)||xf ) ||^ o} J|(a T -^)n||^||T^§o_ i 

+c(e )||^ ||^ o} j|(a r -5 SB )u||l B a||r-||icc (i + M°|||o_ i 

(5.89) 

On the other hand, 

b:= @P ([xf)°]{T T + f x )d x u) - [xt)°]@P ((f T + f x )d x u) 
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and we have 
f x- 2a ~ 1+2 ^b 2 dxdv 

= \\x^9^ ([xt)°]{? T + T x )d x uj - x^xi) ^ ((T r + T x )d x uj \\%>, 

< C||(T r + T x )5 a: n||^a||xb ||L + C'||(T r + f a: )a :E «||| r a \\xi)°\\lo 

k k — 1 {x — 0},1 

< C(e )||T- + T-|||^||a x n||^||^ ||^ + C||xf,°||^_ o}i 

x lll^nll^HT^ + T-ll^ + \\d x uf^_fi T + ^11^^ } 

< C(eo)||T r + T a! ||i „||a a! «||^||xf) ||^o 

U k 

+C( eo )||xf, ||L ||T T + T-||i 00 ^ >A [«(r)] 

+c(||h || L ~)||*h ||L \\d xU \\% s 

{x — 0},l u 

x {|pr + t% Li + ||T- + T»||£. (i + M°|||oJ } . 

(5.90) 

The same holds for the term 

c := 9 P ([xt)°]r A d A u) - [xh ]^ (T a <9 a u) 

and we have 

< C(e )\\T A f LOO \\d A u\\% s \\xi)Y n 

+C(6 )||xh ||L IIT^II^H^nll^ 

{x — 0},1 re — J. 

+C(||f ) ||^)||xh°||2 ||fti«||i»||T A ||L, 

{a;— 0},1 u k—1 

+C( e o)||x^°||^ =o} J|^||^||T A ||| 00 (l + H°II^J • 

(5.91) 

This provides the right estimate for ^4g°, and hence for of A 00 . 

An identical estimate is obtained on the fourth term of the com- 
mutator (5.61). This finishes the estimation of the commutator \O s ,@P]u 
appearing in (5.59), and the proof is complete. □ 
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5.2.4 Conclusion 

The proof of the Proposition 5.7 used essentially Stokes's theorem, the 
weighted Moser type inequalities (A. 34) and (A. 35) of Proposition A. 3 of [13], 
and the weighted substitution inequality type (A. 31) of the same reference. 
One of the points there is that all the constants appearing in these inequal- 
ities are independent of X2 (recall that the sets M X2iX1 there corresponds to 
the sets H\ T here) which is the distance between the boundary of M X2iX1 
and the null hypersurface jV = {x = 0}. So, in our case, all the constants 
involved in the proof of the previous proposition are independent of A. This 
allows us to take the limit as A goes to in (5.52) and obtain an identical 
inequality with E% x [u(t)] there replaced with E%[u(t)]. Therefore we have 
proved the following: 

Proposition 5.10 Proposition 5.7 remains true with A = 0. 

Inequality (5.52) with A = is the key in deriving an existence theorem for 
the Einstein-Maxwell equations with data on a hyperboloid, singular near 
{x = 0}. In this case, we will show that all the and ^ norms appearing 
in this inequality are controlled by the energy. 

It turns out that the proof, in Section 7, of global poly homogeneity of 
the geodesically complete metrics constructed by Loizelet requires a slightly 
different inequality. For this we need to split the metric into two parts as 

Q a(3 = ^ + _ 

The rationale behind such a splitting is, that the Lorentzian metric q will 
be fixed (in fact, it will be the flat Minkowski metric in our applications), 
while the correction 5q will eventually depend on the fields. This leads to 
the obvious corresponding decomposition of T, 

r = f a + 5T a . (5.93) 

We assume that there exist constants a, M and N such that for r € 
[to, ti] we have 

M > ||(^,0«,f)||^ (HT) + ||(^,^,5T)||^ o}i(Hr ) 

+ ||(^-9 r )0 tt ||L-(H r ) , (5.94) 

N > \\(d T u,d x u,d A u)\\gg^ nT ) + \\(3K^)\\Lo-(n T ) 

+ ll(0 rA ,0 xA )ll^_ T(Hr ) + II(^ 5 ^)||^ =o}i(Ht) . (5.95) 

We then have: 
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Proposition 5.11 Let k > n/2+1, a e R, a < —1/2. There exist functions 
C3(eo,Co,a,k,n, M) and C4(eo,Co,a,a,k,n,N), monotonously increasing 
in M and N, which we write as C%{M) and C^{N), such that for all 

T G [t ,ti] 

and for all u satisfying (5.1) we have 

E%[u{t)\ < E%[u(t )}+£ {c 3 (M)(e%[u(s)] + ||F( S )||^« (Ht )) 

+C 4 (iV) (l + imKstfJVWl*-*^) }ds . (5.96) 

Proof: The result is obtained by calculations very similar to those of Propo- 
sition 5.10. We follow that proof until (5.60), which is rewritten as 

□g = °^ V d% + W v d% + + 5T v d v . (5.97) 
This leads to the following rewriting of (5.61): 



--:5A 



■■ M =:ia 



:A 4 



{ & (Sg^dad^u) - 6 S a »$>P (d a d^) } . (5.98) 



:8A 4 



The terms Ai := A\ + 5Ai, i=l,2 are estimated as in (5.62)-(5.63). For A3, 
instead of (5.64) the estimates proceed as before, except that at the end one 
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invokes the weighted Sobolev embedding [13, Proposition A.l]; e.g., 
= \\x^^ (f T (d T - d x )u) - z*TV ((d T - d x )u) ||^ (Ha , t) 



< C s (\\{d T -d x )u\\%^ \\U + \\iPr-W\^U \\U 

\ u re ft— i {x— 0},1 

< cf||f r ||| + ||t T ||2 Wkt)]. (5.99) 

For 5 A3, the following version of the inequalities of [13, Proposition A. 3] 
should be used, for any a,/3,c/ (the proof is identical to that given there): 

\\fg\\^ < C(\\f\\^ 0y j9\\^ + \\9\\^ Jf\\^ + ^) , (5-ioo) 

and, for /3 = (/3i,/3'), 

||x^^(/ 5 )-(x^^/) 5 )||^ + , 

(5.101) 

Instead of (5.64) we then have 



[ x -2«-l+2/?i { S jy2 dx dv 

= \\x^3)P (5T T (8 T - d x )u) - x^ST T 3 p ((d T - d x )u) ||5r «»(H A , T ) 

< C s (\\(d T - d x )u\\%A5T-f + \\(d T - d x )u\\%e« Jl^ll^o ) 

V =*fc k ~ 1 {a;=0},l J 

<c(\\(d T -d x )u\\%«J6r T \\* +\\6T r \\* 9 , ^[n(r)A . (5.102) 



An identical treatment applies to the remaining three displayed equations 
following (5.64). 

The term A$ is split into A^'s as in (5.66), and then for [iv 7^ 00 we split 
j^v _ ^4/^ + in obvious way. All the terms A^ with \iv / 00 are 
then treated as in the proof of Proposition 5.10, and at the end we invoke 
the inequality, for k > n/2 + 1, 

%™ < C\\f\\%>o, . 
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The terms involving 5A^ V with [xv ^ 00 are treated as in (5.102); for exam- 
ple, (5.67) becomes 

[ x- 2a - 1+2 ^{5A AB } 2 dxdu 

= \\X^^ {5Q AB d A d B u) - X^5 9 AB ^ {d A d B u) ||^(H A , T ) 

< C s Y,(\\d A u\\% a -A\8Q%<, + \\d A u\\ 2 ^ k 45Q%, 

A v 1 "J' 1 

< C ll^ll^ll^lll-- + \\58% 0{x=Q} ^ >x [u(r)}) . (5.103) 

In (5.82) it is convenient to use the splitting f) = f) + <5fj. The terms 
involving f) are estimated, using the Sobolev embedding, by E£*[u(t)], while 
for those involving 5t) we write 

= 2||^ (x^°r A ^(5 r - d» - xS^®? {Q TA d A (d T - d x )u) ||^ o a 

< c||g rA a A (a r -^)u||^||x^ ||2 _ CT + c||0 Tj4 ^(a T -a a: )n||| rQ 

^fc fc— 1 {x — 0},1 

(5.104) 

The first line above is estimated as 

C||0 tt |li°°ll(^T -5x)n||^||(50 tt || 2 # ,_ (T , 
as desired. The second is estimated as 

C\\S&% ? [\\{d T - d T )u\\%, W A \\l a ^ + ||(5 r -^||2 llfl"!!^ ) 

{x— 0},1 I 1 fc— 1 K {x— 0},0 J 

^{^=0},! 1 & k-l & {x=0},0 ' J 

(5.105) 

To estimate the term A® (compare (5.83)) we need to split both f)" and g" 
into two. The terms there involving (j" and g^ can be estimated by E^\u(t)\. 
The terms involving 5g^ are estimated as in the analysis of 5 A® . The mixed 
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term involving g" and 5f) is handled in the obvious way 
\\x^@P ([x6l) )[i AB d A d B u]) - x^lxSt) ]®? ([g AB d A d B u] \ \\%ec, 

< C(\\i AB d A d B u\y«\\\\xSb% a -* + \\x8t)°\\<#? i \\i AB d A d B ii\\^ ) 

u k {x — 0},1 re— 1' 

< Cdl^llLoollSAMlbdl^f) !^^-., + ||0 tt |Lo ||dAu||^ Q ||x5f)°|Uo ) 

- 1 - ^fe ft — 1 K {a? — Oj,l 

< c{\\$\\ L ~\\dM\mi¥*% r ° + IIS'II^JI^II^aIuOO]) • (5-106) 
A similar analysis of the remaining terms proves the proposition. □ 

6 Application to the Einstein-Maxwell Equations 
in wave coordinates and Lorenz gauge 

6.1 Change of coordinates 
6.1.1 On the gauge condition 

Throughout this section, the (unphysical) conformally rescaled metric is 
denoted by g, and the (physical) metric is denoted by g; thus = Vfig^. 
Remember that in the original system of coordinates (x^) we have 

UgX^ = with g = rj + h , 

which leads to 

d ll (<r y /\tetgb=0. (6.1) 

We want to rewrite the above equation in the new system of coordinate (y a ) 
(see (6.4)). We have 

V\^9\ = i + \n aP Kp + Q(h) , 

where Q has a uniform zero of order two in h. We set 

In what follows, we use a generic symbol Q for functions which have a 
uniform zero of order two. We have 

fyGTvlditd) = d,[g^{l + ^h a p + Q(h)}] 

= dl ,[{n^ + H^}{l + ^h a p + Q(h)}] 

= d„H^{l + \^Kp + Q{h)} + {rT + H^i^d^Kp + d^Q{h)} . 
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Using this identity, equation (6.1) takes the form: 

(6.3) 

Let us rewrite this equation in the system of coordinates (r, x,v A ) where 

y^= — a g , r = y°<0, x = -y°-p>0 and y l = pJ(v A ) . 

(6.4) 

Recall that 

n = -y a y a = t 2 - p 2 = x(-t + p) >0, (6.5) 

a 71 — 1 

and / = Q 2~ f (not to be confused with division by q tt , as used in the 
previous section), so that 



n^{-(n- i) y , -a^-- iy»y a Q-;}f , (6-6) 



dxv I v dyv % 

thus the left-hand-side of (6.3) can be rewritten as 

We want to analyze the structure of the right-hand side of (6.3). This 
expression is made of three terms which will be labeled Ri, R2, and R3. We 
have (see (6.6) and recall that y a Jp- = 2Q): 

+Q(^k n^y^Hn . (6.7) 
Now, since ^ has a uniform zero of order one, we have 
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Thus R2 reads: 

R2 = ^^^{(n-l^+n^ + Sy/i}^ 



Ti — 1 a n — 1 a n — 1 *> n+1 



=: Q(n^h, Q^y^) + QiQ^H^, n^d^h) + Q{^y^ u , ^y a -^h) . 

(6.8) 

Next 

^3 = -iTd^h) 

= rr |Q(nV^,fiV^) + Q(nV^n^a>) + Q(nV/i,fiV^A^)| .( 6 .9) 

From this, we obtain the following form of the gauge condition (6.3): 

+ fi-V( J R 1 + J R 2 + J R 3 ) . (6.10) 

Now we recall that 

H P> - jf = -f,T + Q^(fr) , 

where h)"* = ffjf^h a p. Therefore 

Equations (6.7)-(6.10) lead finally to the following form of the gauge condi- 
tion (6.3): 

y^-l^(H) = T l_{ !!; |- + 2„^}(^--i^tr,if 

71—1 , 71 — 1 """" Ti — 1 

+n~— Q(n— H, Q—dH) 

+fi-VQ(0¥ff,fi¥j / «_5) . (6.11) 
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We will need the following consequence of this equation: multiplying by y v 
and commuting derivatives one is led to 



(n-5)y u y^ = 2y a A L^H^ + \^ v h\ 



n-1 



Q(n— H,n— h) 

7i — 1 _ , n— 1 Ti+1 „ . 

+n~—Q(n—H, 0,—dH) 

+n^Q{^H,^y a —H) . (6.12) 
6.1.2 On the wave equation 

In wave coordinates (x^), we consider the following wave equation 

In order to check all the hypotheses made on components of the metric in 
our theorem on the energy estimate, we have to rewrite this equation with 
respect the system of coordinates (r, x, v A ) used there. According to our 
previous calculations, equation (6.13) can be written as 

where 

So, let us express the second term of the above equation in terms of coordi- 
nates y v . We already know the identity: 

(6.15) 

with 

^ o (A" 1 = 20^ W + 2fi^7/ A1T y T + 2% tA y a + Srfc^jf j, V 
and 

= tfsffi + 4y W V + 2f](^y A y^ + 5^y a ) . 
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These identities lead to 

H x »V Xfl = H x » fad^y" + 2nS^ e y e + 2n VflX y a + 8 VXaV ^y a y e } |4 

(6.16) 

Now we also know that 

nv/n|4-(n-iWl • (6-i7) 



QyOL Qy 

This implies that 2 

H x »V Xll = 2Q^H X 4 (n-1) {0^ + 2 W } /+(2fW^ A + ft W Q + Ay fl y x y a ) 



df_ 
dy a 
(6.18) 



On the other hand we have 



+ (n - 1) [(n - 3)^^ - ttriap] / j , 
which leads to the following expression of H Xfl K x ^ : 

{d 2 f ( df df \ 

n2 dy^dy~? " (n " l + Va dy~P J + (n " X) [(n " " / 

and after simplifications, we find that 

a 2 / 



dy a dyP 



+ (n - ljn 2 ^^" i 2 (2y W Q + fi^) |^ + [(n - 3) W - Qr^] / 



(6.19) 



2 Note that in this equation, the term y /J .yxH llX is the one which has the the smallest 
multiplicative power of SI. 
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With the expressions (6.18) and (6.19) and writing = Q V# A ^ , equa- 
tion (6.14) reads after simplifications 



d 2 f 

dy a dy> 3 



+ 20^^/ {2(n + l)y M y A y Q + (n + l)fW> + ^Oy } |£ 

+(n - 1) {(n + l)y M y A + fir/ AjU } / 



= ^f (n^/, ^" 1)/2 { - - ~ <» - 1)*}/ 

= : n-^F^/.O^^ . (6.20) 

We want to apply the energy estimates of Section 5.2.3 to the equation 
considered here. So for consistency of notation in that section, we write the 
above equation in the form (recall that = x(p — r)): 



with 



□ git = T(u, <9li) , 

u = f , 



(6.21) 
(6.22) 



{{x(p - t)} 2 5^ x + ^y x y a yP + 2{x(p - r)}(5«y x yP + 5^y a )} , 



(6.23) 



(in order to reduce the typographical length of formulae we will sometimes 
write ip"v for ip al3 ^ u ) and 



du 



n + 3 ~ / n-1 n-1 dlL 

Q-—F Q—u,Q— — 

dy v 



+ |t q - 2n^H x » {2(n + l)y At y A y Q + (n + l)JW£y A + Vx^y a } 
-2(rt - l)n^H^ {(n + l)y^y A + Q^} u . 
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So, we have to check that the metric q defined by (6.23) and the harmonicity 
functions 

TP = J- - d v jvldet^B^) (6.25) 

Vldetgl >■ J 

satisfy the hypotheses of our theorem. 

The tensor ip a ^ ^ defined in (6.23) has the property 

r?a/# a V = (6.26) 
which implies that the contraction 

gains two powers of Q, as compared to a direct power-counting based on 
(6.23). Furthermore, the structure y^y^y^Dv of the term without powers of 
in ip a P ^ v implies that any contraction of the form tp a ^ \iur} a pip po '~g acquires 
an overall multiplicative factor of So if we set 

SQ% ■= Q a ^p ~ ^ , 
it follows that for k > 2 we have 

((5g) fc ) a p := 5 Q a ai 5g a \ 2 ■ ■ ■ = Sl^Q^sf^H) , 

where we use the symbol Q k to denote a smooth function (in this case, a 
polynomial) with a uniform zero of order k, and which may change from line 
to line. A similar analysis shows that, again for k > 2, the trace 

Pk (5g) := ti(5g) k = S S a ai S S a \ 2 ■ ■ ■ S 8 a ^ a = ^l k Q k (^H) (6.27) 

(no summation over k) gains one more power of f2. 
Set 

A a p := 8 a p + 6 S a p . (6.28) 

Equation (6.27) implies 

n 

Pi (A) =tr(J+<5 fl y = J2 C iPj( S 8) = n + l + itiSQ + tfQzin^ H) . (6.29) 

3=0 

Let W(A) denote the characteristic polynomial of A, 

W(X) = det(A - XI) = det A + w x X + . . . + w n X n + (-X) n+1 . 
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Then the coefficients W{ are homogeneous polynomials of order n + 1 — i in 
the entries of A = I + 5q, with w n = (-l) n tvA = (-l) n (n + 1 + tr<5g). It 
is a well known consequence of the Cayley-Hamilton theorem (see, e.g., [26, 
Theorem 1]) that both detA and the w^s can be written as polynomials in 
the p^s, and since each Pi(A) has a factor O 2 in front of the Q2 terms, we 
find that the w^s take the form 

Wi (A) = Wi(I) + £i(tr5 Q ) + D, 2 Q 2 (n^H) , (6.30) 

where £j(tr<5g) is linear in tr<5g. 
Now 

gf = ^ ap ff = (8 a p + 5g a p ) if = A a p if , (6.31) 

hence 

dets 11 = -det(A) , 

which shows that 

detg* = -1 + tt 2 (-Q^tr^H + Q 2 (^^#)) = -1 + ^Q^Q^H) . 

(6.32) 

From the Cayley-Hamilton theorem we have 

A ~ X = ~dkl + " ' + WnAn ^ + (- 1 ) n+lA ") > 
and we conclude that = {rj~ l A~ 1 ) a p takes the form 

Q a p = — — L-n^F^ + n 2 g 2 (o¥ff) 

i + n 2 Qi(n— ) v 

+fi 2 Q 2 (fiV#) , (6.33) 
where the indices on ip a pn V have been lowered with the metric rj a p. 

6.1.3 On the components of the metric 

Recall that, to obtain energy inequalities, our hypotheses on certain com- 
ponents of the metric were 

0°O = _ 1 + xf) O. Q P = _ xt) l. Q 0A + qP A = _ x[) A and gP p = 1 + xf)5 

(6.34) 
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where the functions rj, f)°, \) A are bounded on bounded sets. Since (compare 
(6.4)) 

tf>' = ff*UH, 0^ = ^, 0M = s ii w .^ and 5 pp = ^., 

from (6.23) we have (note that y l Wi = p, p^i^ = y fl + t8V): 

f,° = x ^{p-r)^H^ {{x(p - r)} 2 ^ + 4r 2 W + Ar{x(p - r)}^y A } , 

(6.35) 

f,l = _ x ^( p _ T )^^ _ r)} 2^^ + 4rpy ^ A + 2p{x(p _ T)}yx{5 p + T ^)} 

(6.36) 

f) = x^(p-r)^^{{x(p-r)} 2 ^^o; i a; i +4p 2 y At y A +4{2;(/)-T)}yAP^yA^} , 

(6.37) 

If = _rf( p - T ) V {( p _ r ) + //'■') - 2y A #^} ^ . (6.38) 

We see that the components of the metric (6.23) have the right structure 
(6.34) if the space dimension n is greater then or equal to 7. We will see 
in Section 6.1.1 (see (6.12)) that this can be lowered to n > 6 using the 
harmonic coordinates condition. 
We note the identities, 

a a dv A ^ dv A dv A 

which justify that q oa + g pA has the right structure. In particular, for this 
component the condition n > 4 suffices to fulfill the structure condition. 
We will also need 

g TT = -l + 0(x^), q tx = l + 0(x^), n = O(x^), (6.39) 
s * A = 0{x n r L ), Q AB = V AB + 0(x^) . 



6.1.4 On the harmonicity functions 

Now let us look at the harmonicity functions, defined as 

:= ^— dv {yidet^rj . 
Vl detfl| > 

Since our energy estimates have been established using the coordinate sys- 
tem (x,t, v A ) as defined in (4.10), we need to calculate T M in that coordi- 
nate system. But so far we only have the expression of the metric in the 
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y^-coordinate system. To avoid confusion let us write ^ 2 ^T for T associated 
to the coordinates (r, x,v A ) and ^Y for that associated to the coordinates 
y^. To understand the behaviour of Y under coordinate changes, it is useful 
to write the Christoffel symbols Tp of the metric q in the form 

-pa fa I /~iot 

where the FJg 's are the Christoffel symbols of the Minkowski metric rj, and 
Cg is a tensor. Then, in the coordinate system y M we have 

(i)T Q = -g^ 7 , (6.40) 

=:C a 

since the FJg 's vanish in the -coordinates. Note that C a as defined in 
(6.40) is a vector field, being the contraction of two tensors. In the coordi- 
nates (r, x,v A ) we have 

(2) T a = _ g /3 7 ^ a + Qo \ = a _ Q* _ ( g 41) 



Thus, to calculate ( 2 )Y we need to vector-transform C a to the (r, x,v A ) 
coordinates, and calculate the missing term 0^ 7 f^ above. We start by 
calculating the vector field . We set 

thus 

as in (6.23); we hope that the clash of notation with the completely different 
K Q p appearing in (6.15) will not confuse the reader. 
Prom (6.32) we have (recall that Q means Q2) 



detfl| ) = 1 ± -n—tr v (H) + ST Q(ft— #) 



Thus in the coordinate system y^, 



+n 2 Q^ u (n !1 2 1 H) , (6.43) 
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+n^{(5 - n)y v K^ + nd u K" v } + {o^^fO^ 

and since 

~ y^"" = -0*/^ + 2y a y»} , (6.44) 

and 

= d„H a ^ + 2(n + Z)ypH aP {ft5£ + 2y a y»} + 2fiy"tr IJ if , (6.45) 
we obtain 

^ (g^v/ldetil) = ift 2 ^{(n + 3)i^tr,H-»r'n^t^-ff} 

+ (3n + l)y^(n*£ + 2y a /) 



Multiplying this last identity with yy\ detg|J we then obtain the follow- 
ing expression for the vector field C^: 

+n^{d v H a ^ + (3n + l)y^ {JW£ + 2y a y»} } 
ifi^^oVff) + tfQV^n^H, VL^d u H) . (6.46) 
Now writing the vector field C as 

C = Cd^ =: C T d T + C x d x + C A d A , 

one is led to: 

C T = C° , C T + C x = -LOi(v)C l , C A = — C l . 

dy l 

In order to have all the harmonicity functions in the (r, x, ^^-coordinates, 
it remains to calculate the term 0^ 7 r^ 7 of the formula (6.41). In these 
coordinates the Christoffell's symbol of the Minkowski metric read: 

Kb = 0, 



a/3 

^Tfj, = r x M = 0) ^AB = PXAB 



-pA _ pA _ -pA _ n -pA _ -pA _ xA -pA _ .A 
tt — L tx — L xx — u > 1 tB — 1 xB — ~~°Bi 1 BC — IBC ' 
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where we have denoted the round metric on the sphere by x> an d its cor- 
responding Christoffel symbols 7^ c - These identities lead to the following 
(see identity (6.42)): 

S^f% = (6.47a) 

= pQ AB XAB = !l_± +pn ^H^AB XAB (64?b) 

ifTty, = -- p (Q TA +Q xA ) + 9 BC j£c (6.47c) 
= j 2 C A + + V^c) ; (6-47d) 

where C A = X BC Jbc ^ s m i nus the harmonicity function on the unit sphere. 
Finally, we obtain that the harmonicity functions of the metric g in the 
(r, x, u A )-coordinates read: 

(2) T r = _ c o ( QA8s ^ 

(2) T r + (2) T x = u J v)C i_H^l_ p £l^H^AB (648b) 

P M 
(2)^-4 Qi h(J A 

dy % p 2 

-n^H^ + V^c) • (6.48c) 

We revert now to the notation T for what was denoted by ( 2 )Y above. 
6.1.5 The source term T 

Recall that the source term in ^coordinates reads: 

+ | (1) T a - 2^H X " {2(n + l)y, iyx y a + (n + l)nS^y x + V x^y a } 
-2(n - 1)^^ {(n + 1) W + n^}« ■ 
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From (6.46) we have 

(i) T Q _ 2n^H x » {2(n + l)y^ x y a + (n + l)n<S£j/ A + r, x ^y a } 
= -n 1 ^ | (n - l)y a tr,,# - i] ai/ nd u ti v H^ 

+^{^ x d v H^ + (n - l)y x H» x {Q5% + 2y M y a } } 

This shows that the source term takes the following form: 

_ / \ „ n+3 ~ / n-1 „ n-1 <9li \ 

dy"J V 

^(n-ljUx^KnlljpA + Sl^ii 

{(n - l)y a tr v H - r, av nd u tv v H} |^ 



+0¥{^p^V + ( n - l)y x H^ x + 2^y a } } 



du 
dy" ' 
(6.50) 



6.2 The Einstein-Maxwell case 
6.2.1 Existence of a solution 

The Einstein-Maxwell equations, in harmonic and Lorenz gauge, take the 
form (6.13) (see [4, 16, 18]) with the following replacements there: 

/ = (<V-V,A*) and H^ = g^- V af3 . (6.51) 

Recall that, if v is an arbitrary function, then 

v = n~ n ^ 1 v. (6.52) 

Therefore, we have 

For consistency of notation with Section 5.2.3 we set 

f = u. 
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In this notation 

k k 

and, since 

frap = _ r pm r p V ^ + n -(n-l)/2go0 ^ n (n-l)/2^ ^ ? 

where Q"' 3 has a uniform zero of order two, from [13, Proposition A. 2] we 
obtain that 

< C (j| hf^vW^ooj ||u||^,e-(n-l)/2 

< C Ml^Hiooj . (6.53) 

We define the energy E% x [u(t)] as in Equation (5.47) of Section (5.2.3), 
the metric being defined by (6.23). Recall (see Equation (5.49) of Section 
(5.2.3)) that this quantity controls the i^, Q -norms of df. Now, 

\\dH^\\%e < ||9(7f V"V)II^ + \\9 (^-^g Q/3 (^ (n " 1)/2 V)) \\%,e ■ 
Since 

(n — 1 a , ( i^/n" . \ Ti + 1 n , n — 1 . n— 1 o , n — 1 n— 1 _ , 

J2-— Q a ^(n( n-1 ^ 2 ^)] = Q-—Q a f 3 {Q—h) + n~—Q al3 (n—h,n—dh) 

we have the estimate: 

+l|o- Ii ^ +a Q Q/3 (o Ii ^(/ l ,x- a a/ l ))||^ 

< cdi^n^)^!!^-^,/, 

+c(||&, ar a afc) || L ~ ikii^, s - a ^)|| 9 -„-(„-i )/a 

< C(||A, X-^)|| L oc (||A||^ + ll^ll^-c + ||<%||^) 

< C(\\h, x- a dh)\\ L °o (\\h\\^e- a + HSftllj^) • 

Thus, 

\\dH a0 \\%o < C (\\h, x- a dh)\\ L oo^j (\\u\\^e- a + \\du\\je/) ■ 
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To continue, we suppose that at x = x\ > the maximal globally 
hyperbolic development of the data exists for r G [to,ti], with 

Mi := \\f\ {x=xl }\\L™ < oo . 
We define (compare (5.51)) 

M(r) := ||^||^ (Ht ) + 11(5,(^-^)11^^) + ||(g», h«, Vf^^) 
+ \\{(dr ~ d x )f,dJ,d A f)\\% mr) + \\f(r)\ {x=Xl} \\ L ~ , (6.54) 

with the functions g", f)", T^ 1 = ^ 2 ^T and J 7 defined by equations (6.23), 
(6.35)-(6.38), (6.48) and (6.50). 

For any positive function N(r) we set 

jV(-r) := sup N(s) . (6.55) 

s£[to,t] 

We then have the following: 

PROPOSITION 6.1 Let k G N, a G (—1,-1/2]. Consider the Einstein- 
Maxwell equations (6.13) in space-time dimension l + n>7ifa = — \, 
and 1 + n > 8 otherwise. Let f be defined in (6.51), suppose that to > and 
assume that the initial data, given on the hyperboloid 




(6.56) 



in Minkowski space-time, are such that: 

f |^ o) G {jf k a +1 n (0(^ o )) , and ((^ - dj, d A f) U^ ) € (<M^b)) 

There exists functions C^{n, k, €q, Cq, a, M) and C^(n, k, cq, Cq, a, M ), monotonously 
increasing in M , which we write as C^{M) and C^(AI), such that the energy 
of the system as defined in (5.47), Section 5.2.3 satisfies the inequality 

< 2{M 1 2 + ^[/(r )] 

+ f T C 3 (M(s))E%[f(s)]ds} , (6.58) 

J TO 

where tq = Furthermore, for n + 1 > 7 and a = —1/2 one has 

ll/V)lll~ + £ fc Q [/(T)] < + ^[/(t-o)] + |k^- 7 >/ 2 ^y^(T )||^o 

+ [ T C 4 ( M(s) ) E%[f(s)] ds\ . (6.59) 

Jtq ' 
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Remark 6.2 For n > 7, a prefactor fi - 2~~ in the fourth line 3 of the nonlinear 
term in (6.50) still leads to the estimates here. This remark is important for 
the estimation of the time derivatives in Section 6.2.3 below. 

Proof: For all < x < x\ the trivial identity 

f(T, x) = /(r, xi) - J d x f(r, s)ds 

leads to the estimate (recall that a > — 1) 

||/(r)IU" < M 1+ f 1 \\dj{r)\\ n s«ds 

< M 1 + \\dJ(r)y«. 

From this one easily concludes 

\\f(T)\\yo < C{M 1 + \\(dJ,d A f)(T)y«J . (6.60) 

Now we apply Proposition 5.10 of Section 5.2.3. To obtain (6.58) we will 
show first that, in the Einstein-Maxwell case, the J^, a -norni of the source 
term, the ^-norms of g", and T M are controlled by the energy. Let us 
start with the ^-norm of g*. From the expression of g given by (6.23) and 
the estimate (6.60), if n > 5 then 

||g tt (r)||^o < c(M l + \\dH\\ 2 n ) 

< C (Mi + E% jX [u(t)]) . (6.61) 

The same holds for fj" but with the constraint that the space dimension n 
is larger than or equal to 7. We will return later to the question how to 
improve on the dimension on this term when a = —1/2. 

To estimate the harmonicity functions ( 2 )T given by (6.48), we start by 
estimating the functions C A \ We decompose = + C% + Cg , each 
corresponding to a line in (6.46). The first and second terms are estimated 
as we did for g" and fj": 

||Cf|||o < C^x^H^ + Wx^dH^) 

< C(M 1 + E^ x [u(t)}) for n>3, (6.62) 

3 The fall-off of the component of this term with the lowest power of $1 can be improved 
using the gauge condition. 
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and 

||CXo < C{\\^H\\l„ + \\x^dH\\lo) 

< C{M 1 + E^ x [u{t))) for n>5-2a. (6.63) 

To estimate Cg we recall that its components have a uniform zero of order 
two in H and (H, x~ a dH) respectively, with the second term linear in dH, 
thus we can apply inequality (A. 31) of [13] on the ^-norm with I = 2, /3 = 
We obtain: 

\\cg\\*, < \\Q^n^H)\\l_ 2 + \\Q^(n^(H,x- a dM\\l-?- a 

k s k 

< crcilC^II-f-)II^HSff— C7(||(^, a? -«»^)||i-,)||C^-,^- <,, »^)|lS.»- 

< C (\\H,x- a dH)\\ L oe^ (\\H\\ 2 ^o + \\dH\\l^j for n>3-a. 
Thus we have: 

||CXo ^ C (\\H,x- a dH)\\ L ~) + E^ x [u(r)]) for n > 4 . 

(6.64) 

Note that the function C (\\H, x~ a dH)\\Loo^j will give a contribution to the 

function C2(M(s)) of (6.58). The remaining terms of ^ 2 ^T as given by (6.47) 
are estimated in a similar way. They are controlled by 

C (1 + E^ x [u{t)]) for n > 5 . (6.65) 

We continue by writing the source term T (see (6.50)) as a sum of terms, 
each of the following form 

x"Fi(.,a**(f, x~ a df)) . (6.66) 

Note that all terms are polynomial in df, at most quadratic in df. For 
instance, the first term F arises from products of the Christoffels in the 
Ricci tensor, and from the products of the derivatives dA of the vector 
potential A in the energy-momentum tensor. We then write, for example, 
in the x^ coordinates, 

T 2 ~ (g^dgf ~ F(gt)dgdg = x 2a F{g*){x- a dg){x- a dg) ; 

we then express this in term of h^ u , transform the whole expression to 
the -coordinates, and finally reexpress hn V in term of Huv- This formula 
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Pi 


Qi 


4 


constraint 


n > 


Ti 


n+'i 
2 


^+a 


2 


n > 5 — 2a 


7[6] 


Ti 


n— 1> 


n—b 
2 


2 


n > 5 


6 


Tz 


^ + a 





2 


n > 3 


4 


^4 


^ + 2a 





2 


n > 5 — 2a 


7[6] 


^5 


2-^ + a 


n— 5 
2 


3 


n > 3 


4 


2-^ + 2a 


n— 1> 
2 


3 


n > 3 — a 


4 



Table 6.1: Restrictions on the dimension from the source terms. 

shows that the T 2 in the Einstein equations have a uniform zero of order two 
in (f,x~ a f). A similar analysis applies to the contribution of the Maxwell 
fields to the Einstein-Maxwell equations. 

We use the following estimate to show that the J£jL a -norm of J- is con- 
trolled by the energy of the system: Suppose that T% has a uniform zero of 
order li in (it, x~ a du), then applying to this function the second part of 
lemma A. 2 of [13], for 

Pi + Uqi > a . (6.67) 
We choose e > so that pi + tiqi > a + e, and write 

||ajWJi ( . ,x qi {u, x~ a du)) 
= \\Fi{. : x^{u,x- a du)) u 



2 

a- Pi 
k 



< C(\\(u, x- a du)\\ L ~)\\(u, x- a du)f K _ v ^ m 

< C(\\(u, x- a du)\\ L oo) (||it||^- £ + ||x- a 9n||^_ e ) 



< 



C(||(u, x- a du)\\ L ~) +£E>(t)]) . 



The analysis of the nonlinear terms (6.50) along those lines gives the follow- 
ing table: Here the -7-i's, % = 1, ... ,4, correspond to the i-th line of (6.50), 
while the two rows for J-$ correspond to the two respective terms in the last 
line of (6.50). In the last column the number in square bracket is obtained 
by estimating below the non-linearity in a more efficient way. 

It turns out that the threshold on the space dimension n can be lowered 
to n = 6 for the components T\ and J-4 of the source term T . The quadratic 
terms in those expressions with the lowest powers of fi are of the form 
Q !1 2~G(rL :L 2~ f)dfdf for F\ and tl^Hdu and Q^dHdu for One can 
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n — 5 

estimate the J^-norm of fi - ^ - Hdu using instead (5.100): 



n — 5 

2 



< C\\\H\\% \\du\\ 2 _ „_ 6 +||xVH||2 iii 



< C(||n||l oo + ||9n|| 2 ^)(||u||^o + ||0n||^) if n>5 



< , c(|Hiioc + ||a«||^)(i + 



V 11 " -^h 



see (6.60) 

< CdMHoo + H&tH^ja) (l + llauH^a) , (6.68) 



for k >n/2. Next, 



A4» 



< c ( ||a#||| a ||an|| 2 _„_ s +||a#" 2 



" " ^ 



\.m< 



< C\\du\\%« Wduf^c, if - ^— — - a < i.e. n > 5 - 2a 

< C||a#||^.E? ifc [u(r)] , 

and so the last inequality will be true provided that 

n > 6 if a = —\ 
n>7 if - 1< a < -\ 

A similar calculation applies to T\. 

These estimates and the table show that 

\\F(u,du)\\^<C(\\u\\ L ~A\Mh «){l + E^ x [u(T)}) (6.69) 

for 

J n > 6 if a = — | 

\ n > 7 if - 1< a < -| 

Inserting inequalities (6.61)-(6.63) and (6.64)-(6.69) in (5.52) of Section 5.2.3 



gives (6.58). 

Now, at several places of the calculations above the term 

V> := y a y P H aP 
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is the one that occurs with the lowest power of f2. It follows from the wave- 
coordinates conditions that this term solves equation (6.12), which can be 
written in the form 

-y a d a ^ + = C , (6.70) 

where 

C := 0(^^(5) +Ife A(^ + ^tr,(H))) 

n — 1 „ , n— 1 7i — l . 

„ — 1 , u — 1 „ n+1 „ ^ 

Q(Q— if, n—dH) 

n — l , n — 1 ri — l „ (3 - — ■ 

+n~— Q(n— h, n— y a —H) 

=■■ C1 + C2 + C3 + C4, (6.71) 

where Q corresponds to the i-the line. The point is that all terms in £ 
contain effectively multiplicative powers of 0. 

Solutions of (6.70) take the form, for tq < r < t\ < 0, 

^(r, x) = (-r)-(- 5 )/ 2 (£(- S ) (n " 7)/2 C (*, f ) + (-ro) ( "" 5)/2 ^ (r , ^) 

(6.72) 

This gives immediately, for any 7, 

-)ll^< ||V(to)||^ + C(^o,ti) [ T \\«s)\\^ds , (6.73) 

J TO 

similarly for or ^ 7 -norms. In the notation of (6.55) one thus finds 

U{r)U! < \mr )\\^+C(r ,n) f T \\as)\\^ds 



TO 



< \Mt )\\^ +C7(^,t 1 )(tl -7a) ||C(r)|| y 7 ■ 

Using this to estimate fj° we obtain 

< C( ||x("" 7 )/V(r)||^o + ||x("- 5 )/ 2 ff(r)||^o 

< C(||x("- 7 )/V(r )||^ + ||x("- 7 )/ 2 C(r)||^o + |M"- 5 >/ 2 £(t)||^o ) 
We have, for example, 

lk (n - 7)/2 Ci(T)||^o < C (\\x^ 2 H(r)\\^o + \\x^/ 2 dH(r)\yo) , 
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which, for n — 5 > —2a, can be controlled by ||.ff (t)||l°° and E^[u(t)] in 
view of (6.60). This requires n > 6 if a = —1/2, or n > 7 if a G (—1, —1/2). 
An estimation of the remaining £j's along the lines of those already done 
above presents no difficulties. 

The functions h 1 and f) have the same structure and so the same estimate 
applies; the function i) A has a higher multiplicative power of Q so that the 
original straightforward estimate applies. 

The final inequality (6.59) follows immediately from this and from an 
obvious version of the estimate (6.58) for the remaining terms in the equa- 
tion. 

We finish this proof by noting that the above treatment of y a ypH a ^ can 
be used to improve the threshold on dimension for some of the entries of 
Table 6.1; this will, however, not improve the threshold on n of the theorem. 

□ 

We are now ready to prove existence of solutions in weighted Sobolev 
spaces. For s > consider the family of hyperboloids: 

*. = {W):J-s = V* + W}- ( 6 - 74 ) 
Let (ft be defined in (3.2). We have the following 

Theorem 6.3 (Propagation of weighted Sobolev regularity) Suppose that k > 
[^j + 1, with n = 6 and a = —1/2, or n > 7 with a G (—1,1/2], and let 
to > 0. Suppose that 

f\ H ^ o) G ( 4 a + i ni °°) (<K^o)) , (d T f,dJ,d A f) |^ o) G JW(^o)) 

(6.75) 

where f and f are defined by (6.51)-(6.52). In the case a = —1/2 andn = 6 
assume moreover that 

.-VV^I^G^ . (6.76) 

Then there exists i* > to an d a solution of (6.13) defined on U S? s such 

se[to, t*] 

that, Vr G [— 2^,— ^7] =: bo) r *] we have: 

f G L°°([t ,t*], JC(H t ) nL°°(H r )) , (6.77) 

(d T f, dj, d A f) G L°°([to,t*], ^T(H r )) . (6.78) 

Moreover, any solution for which M(t), as defined in (6.54), ^ s bounded on 
[tq,ti] satisfies (6.77)-(6.78) with r* =t\. 
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Remark 6.4 Using the weighted Sobolev embedding theorem we conclude 

/V)e (^_[f]-i nL °°)( H -)' ( 6 - 79 ) 



d T f(r),d x f(r),d A f(r)) G ^_ [§] l (H T ) . (6.80) 
when the prescribed data are as in Theorem 6.3. 

PROOF: In order to apply the Gronwall-type Lemma 5.2 of [13], we need to 
prove that all the norms in M (see (6.54) and (6.58)) are controlled by the 
energy or the L°°-norm of u. Since k > [^] +1, from the weighted Sobolev's 
inequality, we have: 

\\(d T -d x ,d x ,d A )f\\%° < \\{d T -d x ,d x ,d A )f\\\ a <El x [u{r)\ . (6.81) 

Let us look at the L°°-norm of (d T — d x )Q^- Recall that the expression of 
g" is given by (6.23). We estimate here only its worse term which is of the 
form Q^2~ H. We have: 



\\(d T -d x )(n-^H)\\ 2 Loo < c(\\n—(d T -d x )H\\ 2 L oo + \\n—H\\l 

< C(\\u\\loo+E^ x [u(r)]) for n>7. 

Thus, 

\\(d T - d x )Q*\\ < C (\\u\\ 2 L oo + E% x [u(t)}) for n>7. (6.82) 
Next, since k > n/2 + 1, as in (6.61), we have 

llJtl|2 s II„I1||2 

W\\^ < WWcq 

< C (Mi + E% x [u(t)]) , for n > 5 . (6.83) 

Similarly, 

\\t)%o < \\i) S \\lo < C (Mt + E% x [u(t)}) for n>7. (6.84) 

1 k 

If a = —1/2 the threshold n = 7 in (6.82) and (6.84) can be lowered to 
n = 6 by using the estimate (6.73) on the slowest decaying term ip. 

To estimate the ^"-norms of the harmonicity functions, we use again as 
in the previous estimate the Sobolev inequality and obtain a control of these 
norms by the energy with the same constrains as in (6.62)-(6.65). Let us 
estimate now the L°°-norm of u. Integrating backward along the integral 
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curve of the vector field Y v d u = d T — d x we can write the identity (here we 
omit the variable v A ) 

u(t, x) — u(t , t — To + x) = I (d T — d x ) u(s,t — s + x)ds . (6.85) 

J TO 

Thus we have 

\u(t,x)\ < \u(tq,t - to + x)\ + / | (t - s + x)~ a (d T - d x ) u(s, t - s + x)\ (r - s + x) a ds 

J TO 

< |u(t ,t - t + x)\ + / || (d T -d x )u{s)\\<#* (t- s + x) a ds 

J T 

< \\u(t )\\l°°+ f ' \\{d T -d x )u(s)\\^{T-s) a ds . 

J TO 

Since k > | we can now write (—1 < a < —1/2): 

IKt)||l°° < \\u(t q )\\l°° + \\{d T - d x )u(s)\\^ (t - s) a ds 

J T 

< \Ht q )\\ L oo+I Je* x [u(s)] (T-s) a ds (6.86) 

Inequalities (6.81)-(6.86) show that from (6.58) we have the following: 

IKt)|||~ + E% jX [u{t)] < C {\\ufo)\\ L - + ^, A [n(T )]) 

[ T $ (EZ x [u(s)\, \\ u (s)\\ L ~)(l + (T-sr)ds, (6.87) 

J TO 



+ 



where $ is bounded on bounded sets. We can now apply Lemma 5.2 of [13] 
and obtain that there exists a time To < t* < depending on ||u(to)||l°° + 
E% x [u(to)] and on the function $ such that Vt G [to, t*] , 

||u(r)|Ucc +E% x [u(t)] <1 + C{\\u(to)\\l~+E% x [u(to)]) , (6.88) 

which provides the desired bounds. 

If one knows a priori that M(t) is bounded, (6.87) becomes effectively 
a linear inequality, and the claimed global bound immediately follows. 

Actually, the solution constructed here is defined on (see Figure 6.1). 
In order to obtain a solution in a whole neighborhood of the hyperboloid JT'o, 
we proceed as follows: Let R > be a real positive number such that the 
level set r = R lies in the region where the energy estimates above apply. 
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We consider the Cauchy problem for (6.13) with initial data obtained by 
restriction on 

y (R) = y n {(^) -.o<\x\<r} . 

We thus obtain a Cauchy problem on a compact region. We can now apply 
to this problem the conclusion of Proposition 3.2, p. 378 of [24]: there exists 
a time t + e]to, 0[ and a smooth solution on (see Figure 6.1) 

y + = u <K^(i?))n^+(<K^o(#))) , 

te[to-2^] 

where @ + denotes the domain of dependence, and where 

S%(R) =S%n {(a^ 4 ) : < \x\ < R} . 
From uniqueness in Proposition 3.2, p. 378 of [24], we conclude that the 




Figure 6.1: The sets V + and % 



solutions constructed on Y + and coincide on^+ n "W Tt which is not empty 
for R large enough. We thus obtain a solution of (6.13) with (6.51) in a whole 
neighborhood of S^q. □ 

6.2.2 Space-regularity of the solution 

For smooth initial data the solution constructed in the previous section is in 
C°°(y+ U % Tt ). In this section we want to show that, for data given in the 

space n we can control the growth, near x = 0, of all space derivatives 

fceN k 

of the corresponding solution. We have the following: 
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Theorem 6.5 Under the hypotheses of Theorem 6.3, suppose moreover that 
the initial data given on the hyperboloid S^q satisfy 

f U^ ) e (^C n L°°) (H ro ) and df |^ o) G ^C(H To ) . (6.89) 

If a = —1/2 and re = 6 we afeo suppose that (6.76) holds for all k. Let r* 6e 
as in Theorem 6.3 with k = ko, where ko is the smallest integer larger than 
[n/2] + 1. Then 

Vre[r 0) Tj /(r) G n L°°)(H r ) , 9/(r) e JC(H T ) . (6.90) 

Furthermore, any solution with smooth initial data as above for which M{t), 
as defined in (6.54), * s bounded on [to,ti] satisfies (6.90) with t* = T\. 

Proof: We provide the details for n > 6; the treatment of the case re = 6 
is similar. From Theorem 6.3 there exists a time r* and a constant C* 
depending on ko such that Vr G [to, t*[ , 

Kr)||i»+^>(r)]<C7*. (6.91) 

Now let /c G N, k > k , since / | 0( ^ o) G (^ Q n L°°) (H To ) inequality (6.58) 
holds. Now the function C3(M(s)) appearing in this inequality is controlled 
by E% \[u(t)] and thus by C*, therefore, from (6.91) we have: 



EI x [u{t)]<C(C*) (l+ [ T El 



[u(s)]ds 



Applying Gronwall's inequality we obtain: 

EtM^ ^ Ce ° Tt ■ 

This inequality shows that, for all k, 

dueJt? k a , (6.92) 
as desired. □ 

6.2.3 Estimates on time derivatives of the solution 

In order to estimate the time derivatives of the solution, we introduce a new 
set of variables (y, x) (compare Figure 6.2): 



T 

X 



? + which implies that 

X I O x — Ox o O r 
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(x = 0,t = t*) (x = cr(r*), r = r*) 




(x = O l2 / = 0) 

Figure 6.2: The variables (x, r) and (x,y), with T := r* — To- The function a 
has been introduced in (5.17). We hope that the reader will not get confused 
by the fact that the boundary x = 0, at the left-hand sides of the figures 
here, is depicted at the right-hand side of Figure 5.1. 

Note that in these new coordinates, the hyperboloid S^q is represented by 
the set {y = x}. Since we are interested in the behavior of solution in a 
neighborhood of the set {x = 0}, as in [12] we restrict our attention on the 
subset of defined by: 

W = {(y,x,v A ) : < x<y, v 6 ,0 < y < 2(r* -r )} . 
Recall that the definitions of the spaces 

^{x=0},fc(^0' < ^{y=0},fc('^')' ^{0<x<y},k{^\ and ^{0<x<y},k(ty ' 

can be found in Appendix A.l page 116 of [12] with d x there corresponding 
to dx here. 

Remark 6.6 In the coordinates (y, x) the components of the inverse of the 
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metric read (compare 6.39): 



g w = 4 (fT + tT) + s xx = o(x 



n — b . 
2 



■iyA 



2q xt + q xx 
2q tA + Q xA 



0{x 



n-l . 
2 



(6.93) 
(6.94) 
(6.95) 
(6.96) 

0~" • (6-97) 

Recall that the hypersurfaces 5? s have been defined in (6.74). As a first step 
towards proving propagation of poly homogeneity, we obtain some informa- 
tion about the ^-derivatives of the fields: 

Theorem 6.7 Suppose that k > [^] + 1. Under the hypotheses of Proposi- 
tion 6.1, there exists t* > to an d a, solution of (6.13) defined on U 



such that: 



f G ( ^{0<x<y},k- 



d T f,d x f,d A f)G^ ^ y} ^ hl m 

where f and f are defined by (6.51)-(6.52). 



se[t , u] 
(6.98) 
(6.99) 



Proof: The proof of existence is given by Theorem 6.3 and we have / G 
L°°(ty), df e ^ =0}fe _[zii We note that from ( 6 - 4 ) and ( 6 - 5 ) we 



have 

n 



-V - 2r ), ydyVL 



-xy, 



and 



d A n = . 

(6.100) 

Identities (6.100) show that if we apply to (6.21) the operator (d A ,xd x ,yd y ), 
then we obtain a wave equation with (u, d A u, yd y u, xd x u ) as the new un- 
known functions in which the coefficients have the same powers of x as in the 
original equation, and the source term the same structure. More precisely, 
set 



U 



( u \ 

d A u 
xd x u 
\ydyuj 



we thus obtain 



U 
dU 



( u \ 

d A u 

xd x u 
ydyu 
du 
d(d A u) 
d(xd x u 
\d {yd y u)J 



(6.101) 
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and let us derive a wave equation on U. Straightforward calculations lead 
to the following identity (here we write the source term as a function of 
variables p\ and pf, that is F = J-~(-, pi,p% ) ): 

n s (yd y u) = -{yd y ^)dl p u + 2Q^d a d y u - (yd y T a )d a u + T^dyU 

+(yd y J r )(-, u, du) + (yd y u) — (-, u, du) 
dJF 

+ (d y (yd a u) - 6'^dyu) —(-,u,du) . (6.102) 

We write 

{ydyz w )%u = d y Q yy (dy(yd y u) - d y u) ~ n^(du + u)du , 

{ydyQ^dydiU = dyQ VV dy(xd~ X U) ~ 8U dU , 

{yd y Q™)diu = 0{x^) {d s (xd £ u) - d s u) ~ tl^UdU see (6.39), 



gw^« = O(xV)! (dy(yd y u) - d y u) ~ Vt^UdU , 

and 

2 B ay d a d y u = flwaju-l^flfu + 2 £A ^a A n + g AB d A d B u + T a d a u - T(u, du)} 

All the terms arising above have a structure similar to (6.50). A similar 
comparison of the remaining terms shows that we have 

D g (yd y u)=F 1 {U,dU) , (6.103) 

where the source term T\ is of the general form as in (6.50) with the differ- 

n — 7 re — 7 

ence that it has a term SZ 2 {Jot/ with a multiplicative il 2 ; this term can 
be estimated as in (6.68) as long as n > 7. Moreover, it is easily checked 
that this remains compatible with the estimate of Proposition 6.1 (see Re- 
mark 6.2). Note that the procedure above introduces into the coefficients 
of the source terms the function (y,x) 1 — > ~, which is bounded on % ; fur- 
thermore, xd$- = —ydy^j = -, which implies that we will not loose the 
regularity of the source terms, as needed for the problem at hand, when 
iterating the process. 
From the identities, 

n s (xd iU ) = -{xd^dlpu + 2g a£ d a d iU - (xd £ r a )cu + 

dT 

+(xdxJ c )[-,u,du) + dr J .(xdxu)——(-,u,du) 

dpi 

QT 

+ (dz(xd a u) - SZdiu) —(;u,du) , (6.104) 



SO 



O g (d A u) = -{d A Q af} )d 2 al3 u- (d A T a )d a u+d A u^(u, du)+dd A u^(u, du) , 

dpi dp 2 



(6.105) 

we deduce that the same analysis holds for U g (d A u) and O s (xdxu). There- 
fore we have derived for the new unknown function U a wave equation of 
the form (6.21), i.e.: 

U S U = $(U, dU) . (6.106) 

In order to apply to this equation Theorem 6.3, we have to check that the 
initial data for U are in the right spaces. Note that the initial data are 
prescribed on the subset {x = y} of . We denote this hypersurface by So, 
thus Sq = <ft(^o) n , and we set 



s s = <t>(y s ) n % c h_ 1/2s . (6.107) 

We want to prove the following. 

Lemma 6.8 Under the hypotheses of Proposition 6.1 we have: 

(u,d A u,xdzu,yd y u)\ Xo G (JC n L°°) (S ) , (6.108) 

(du,dd A u,d(xd iU ),d(yd y u))\z g ^i(s ) • (6.109) 

Proof: By assumption, we have 

u| So G (^ n L°°) (S ), and (8 A u, 8 & u , d y u)\^ G ^T(S ) . (6.110) 
Now, using Sobolev's embedding theorem, we have 

x- a (d A u,d s u ,dyu)\z o G L°°(S ) . (6.111) 
This leads to the following estimates: 

\xd x u\ So \ = x 1+a \x~ a dxu\x \ < oo, 

\ydyU\-E \ = \xd y u\T, \ = X 1+a \x~ a dyU\z \ < OO . 

To see that d A u(tq) is in L°°^q), we proceed as follows: integrating d A u(ro) 
in x until xo gives the inequality 

l-XO 

d A u(r , x ,v A ) - d A u(T , x, v A ) = \ d x d A u(r , s, v A )ds , 



s a ds 



which leads to the estimate 

\8 a u(t ,x,v a )\ < \d A u(r ,x ,v A )\ + \\ d xu(ro)\\^ £=ohl J 



< \d A u(r ,x ,v )| + \\ds;u(T )\\je>c* / s a ds 



T 



81 



(recall k — 1 > §). Since 

\\d A u(T ,xo,v A )\\ L oo^ < oo, ||<9 £ u(t )||^ < £^[ii(T )] < oo , 

and 

f x ° 1 

/ s a ds = (a£ +1 - x a+1 ) < oo , 

we conclude that HSa^C 7 !))!!-/: 00 < °°- Thus (8au) |s £ ^°°(So) an d we then 
obtain (6.108). On the other hand we have 

||(? v (i9xu)|s || i ^_ i ( Eo ) < ll^^^^lsoll^r^^So) + H^^ in l s o||^_ i (s ) 

< ll^ n |soll^(S ) 

< oo see (6.110) . 

Similarly, we have d{yd y u)\-£ Q -, ddAu\% S J i ^_ 1 (T i q). We thus obtain 
(6.108) and the proof of the lemma is complete. □ 

Now, we apply Theorem 6.3 to (6.106) and obtain that 

(u, Oau, xd x u, ydyu) € L°°([0, 2(r* — to)], {J? k a n L°°) (S £ )) , (6.112) 



(du,ddAU,d(xd £ u),d(yd y u)) e L°°(J0,2(t* -T )], ^ Q (S £ )J . (6.113) 

Using once more the Sobolev embedding theorem, we obtain that Ve £ 
[0,2(t* -t )] 

(u,d A u,xd x u,yd y u)\^ e ^^^^(E^ . 
||(u,5 j4 n,x5 £ u,yaj / n)|| L( x )( ^ ) = sup || (u, d A u, xd^u, yd y u)\^ T \\l°°(3? t ) 

r6[ro,r*] 

< sup || (u^A^xd^ydyu)]^ \\^{^ T ) 

t£[tq,t,] 
OO . 

see (6.112) 

Using now (6.113) instead of (6.112) we have 

\\(u,d A du,xdxdu,ydydu)\\ L <x>(w) < oo . 

This allows us to conclude that (u, du) is in ^m <£<2y | \{^)- Now, if we 
repeat this process j times with j = k — [^] — 1 then we obtain that u is in 
^{o<£<2/} fc-2-i(^0- This completes the proof of Theorem 6.7. □ 
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COROLLARY 6.9 Under the hypotheses of Theorem 6.5 we have the follow- 
ing: 

f G K Q o<z<y},oo W n W and df G ^f <*<, } ,oo W • 
Proof: The result is a combination of Theorems 6.5 and 6.7. □ 

7 Polyhomogeneous solutions of the Einstein-Maxwell 
equations 

Let 5 be a positive real number. We recall that the spaces of polyhomoge- 



neous functions £/{ x=0 y, £?f x=0 \, ^{o< x <y} an( i ^{o< x <y} are defined in [12, 
Equations (A.1)-(A.2)]. We consider the Cauchy problem for the Einstein- 
Maxwell equations (6.13) with (6.51) in wave coordinates (x^) and Lorenz 
gauge with prescribed data on the hyperboloid J5*o (see (6.56)) at the in- 
terior of the future light-cone with vertex the origin of coordinates. The 
coordinate x in which the polyhomogeneous expansion is taken is x = 

n 

where t = x° and r = \x\ = ^(x 1 ) 2 ). Indeed we have (see (6.4)): 

i=l 



-T - p 



j2 ( x *) 2 



t r 



-t 2 + r 2 t 2 
1 



t + r 

We want to prove that, polyhomogeneous initial data for the above Cauchy 
problem lead to polyhomogeneous solution. We have the following: 

Theorem 7.1 Consider the Einstein-Maxwell equations on R 1+n , n > 8. 
Let S G R be such that 1/(25) G N when n is even and 1/5 G N when n is 
odd. Suppose that the initial data for (6.13) in wave coordinates and Lorenz 
gauge are polyhomogeneous on the hyperboloid S^q : 

fU G * V <=o } n d T f\^ o g ^< =0} , (7.1) 

with f = {g^u — rj^ u , ) . There exists a time t + > to and a solution defined 

on U S^t such that Vi G [ioi*+] we have: 

te[to,t+] 

f(t) = f\ G x^^j and d T f{t) = d T f\ G x^" 1 ^* . 

(7.2) 
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Moreover, the solution is polyhomogeneous at J? , in the above polyhomo- 
geneity class, as long as it remains in J^(H. T ), for some a G (—1,-1/2]. 



Proof: Choose any a < 0; we then have the inclusion £?{ x=0 -y(<t>(<?b)) C 
^C(0(^o))- It follows from (7.1) that we have: 

/|^ o) G(^CnL°°)(^o)) and df |^ o) G ^C(0(^o)) ■ (7.3) 

For definiteness set a = —1/2. From Theorem 6.5, there exists a time 
and a smooth solution / of (6.13)-(6.51)-(7.3) defined on such that 
Vr G [tq,tJ, /(t) G ^"(H,-). Next, applying Corollary 6.9 one obtains that 



/ G {^{0<,<vhoo n i°°J (*) and 9/ G «?o<^v},oo(*) ■ 
From Theorem 3.1 of Section 3, with 

i>\ = f , V>2 = (0 y f, d A f) , f = d x f , 

we obtain (7.2), and the proof is completed. □ 

It is natural to find conditions which guarantee that solutions remain in 
weighted Sobolev spaces on hyperboloids, and hence remain polyhomoge- 
neous if the initial data are. One such criterion is provided by the following: 

Theorem 7.2 Suppose that k > [§] +1, with n = 6 and a = -1/2, or n > 7 
with a G (—1,1/2]. Solutions of the Einstein-Maxwell equations remain in 
J^? 1 , a G ( — 1, —1/2] as long as f remains in ^r^—gi 1; with 

« > ■ P.4) 

The same is true for 

k > — 2 ^ provided that \\x 2 y ll y v H IJ ' v {tq)^^ < 00 . (7-5) 

In particular, in dimensions n + 1 > 9 i/ie small data solutions of [18, 19] 
evolving out from data stationary outside of a compact set are polyhomoge- 
neous. 

Proof: We want to use Proposition 5.11 to show that solutions as above 
remain in a G (—1,-1/2]. For this, consider first the right-hand side of 
(5.94). For k > — (n — 5)/2 one immediately finds that ||<5fl"||<g'0 is finite, 

{x= 0},1 
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similarly for (d x — d T )5<$ when n > — (n — 7)/2. Finiteness of ||<5fj"||~>o 

{a:=0},l 

is straightforward for re > — (n — 7)/2 from (6.35)-(6.38). The estimate on 
<5Y follows from (6.46) and (6.48) provided again that k > (n — 7)/2. 

For re > — (n — 5)/2 the slowest decaying terms in f), T, and in (<9 X — 5 r )g" 
are handled by the < ^° a . =0 | 1 -spaces equivalent of (6.73), 

{n-7) 

0} ,1 

(n-7) /" r (n-7) 

< |x 2 ^( To )|| o + C(r ,T 1 )/ C(s)|«y ds , (7.6) 

{a;=0},l y {a; = 0},l 

(n-7) 

under the supplementary condition that \\x 2 ^(ro)||<^o is finite. 

{x— 0},1 

For any cr such that 

ct < k (7.7) 

we have 

Hence the right-hand side of (5.95) is finite for all such cr's, and so (5.96) 
applies. It remains to show that the integrand in the second line of (5.96) 
can be bounded by a multiple of the energy: 

\\(5 Q KS^,ST)\\l r , {HT) < CE^u(s)} . 

This is easily checked to hold under (7.4) or (7.5) if we choose a so that 

n — 7 

This, together with (7.7), explains (7.4). 

The property that the solutions of the Einstein-Maxwell equations con- 
structed by Loizelet are in ^ x=0 ^ 1 on all hyperboloidal slices has been 
verified in (4.22). There -re = 5 E (0, 1/4) □ 
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